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(54) Channel estimation using a probability mapping tabic 



(57) The invention relates to an estimation appara- 
tus of a receiving device for estimating the transmission 
operation of a transmission channel, in which case the 
transmission channel is subject to interference, and In 
which the analogue data (5) received from the transmis- 
sion channel is assigned an estimator sequence (a) of 



digital data values, in which the estimator apparatus pro- 
duces at least one estimation sequence (a), In which 
both characteristics of the transmission channel and 
characteristics of the receiver are Included, in order in 
this way to improve the reception characteristics in gen- 
eral and, in particular, the bit error rate of the estimation 
sequence (a). 
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Description 

[0001] The invention relates to an estimation apparatus for estimating the transmission operation of a transmission 
channel for digital data, in which case the transmission channel is or may be subject to interference, and to a receiver 

s and a system for transmitting data using such an estimation apparatus. 

[0002] Optical long-distance communications systems with transmission rates of 10 Gbit/s or more and which use 
standard monomode optical fibres, erbium-doped fibre amplifiers and noncoherent optical receivers, are adversely 
affected by nonlinear pulse dispersion and signal-dependent non-Gaussian noise. If receiver circuits corresponding to 
the prior art are used here, losses occur with regard to (he range, the data rate or the error rate. This is caused by the 

to square-function characteristic of the optical/electrical transducer and a number of physical phenomena such as chro- 
matic dispersion, polarization-mode dispersion, the nonlinear Kerr effect, chirp, extinction ratio, self-phase modulation 
or noise resulting from spontaneous emission, see, for example, Agrawal, G.P. in Fiber-Optic Communication Systems. 
Second Edition, New York: John Wiley & Sons, Inc., 1997. 

[0003] The wide range of effects, some of which vary with time, leads on the one hand to the impossibility of producing 
is a closed a priori description of the behaviour of the resultant electrical channel and of the resultant noise. On the other, 
this results in a nonlinear transmission channel which varies with time and is subject to memory and which has signal- 
dependent noise with an unknown likelihood density, which must be estimated while data transmission is taking place. 
[0004] In addition, problems occur in the provision of suitable receiver circuits for data rates of 10 Gbit/s or more. 
Furthermore, some digital devices operate in such a way that an analog/digital transducer (ADC) is required for these 
20 sampling frequencies. Such components are not prior art and at the moment can be produced only if very poor resolution 
is accepted, with further static and dynamic nonilinearities, that is to say nontinearities which are dependent on previous 
samples. 

[0005] Since virtually distortion-free transmission can be assumed at data rates of less than 2.5 Gbit/s in optical long- 
distance systems, receiver circuits comprising a low-pass filter, automatic gain control and a threshold-value decision 

25 maker are regarded as prior art, by means of which the required bit error rate of less than 10 12 is achieved. If the data 
rate on these optical paths is increased, then, however, such arrangements do not take account of the pulse dispersion 
effects described initially and which now occur. The error rate can then well exceed the permissible limit, depending 
on a large number of influences, as a result of which the connection will fail at times or entirely. Alternatively, error- 
correcting coding would have to be provided, with the usable data rate thus being reduced or the distance being short- 

30 ened. 

[0008] Furthermore, approaches are known for compensating for some of the effects by means of optical equalizers, 
see Helsmann, F. et al. in Automatic Compensation of First-Order Polarization Mode Dispersion in a 10 Gb/s Trans- 
mission System. Proc. of ECOC98, pages 529-530, Madrid, 1 998. A first disadvantage is that the active optica! com- 
ponents required to do this represent a considerable cost factor and, secondly, the signal-theory causes of the problems 
as resulting from this can be only partially solved, and only in comparatively good conditions. The latter also applies to 
the adaptive electrical filter approach (linear feed forward equalizer). This Is because, firstly, this generally increases 
the noise level and, secondly, it can remove distortion only from linear channels or, in some cases, slightly nonlinear 
channels. 

[0007] The structure of the decision feedback amplifier which is known from Kasturia, S., Winters, J.H. in Techniques 
<*0 for High-Speed Implementation of Nonlinear Cancellation. IEEE Journal on Selected Areas in Communications, Vol. 
9, pages 711-717, 1991, is also used for this purpose. This interference compensation principle is, however, funda- 
mentally lossy, and leads to error propagation. Furthermore, the achievable error rate is always dependent on the 
decision time within the bit interval. 

[0008] The theory behind the receiver for nonlinear channels with memory, which makes decisions on the basis of 
is minimum error likelihood, is discussed in the literature, see Benedetto, S., Biglieri E., Castellani, V. in Digital Trans- 
mission Theory, New Jersey, Prentice-Hall, 1987. This solution approach relates to a probabilistic estimator for the 
transmitted bit sequence, which requires a priori knowledge of the signal forms of the nonlinear channel not subject to 
interference, as well as the conditional likelihoods of the interference. The above statements have shown, however, 
that neither is known. 

so [0009] With regard to analogue/digital transducers, known arrangements already exist which comprise a nonlinear 
Instantaneous function, the ADC and a digital correction apparatus, see DE 44 09 063 C1. Since these arrangements 
are used only to compensate for static effects (for example to Increase the drive range), the reaction to dynamic changes 
of channel characteristics or interference Is extremely unsatisfactory. Furthermore, such systems require a fast ADC 
which, as a result of the problems explained, leads to a further adverse effect on the resultant transmission channel. 

55 [001 0] The invention is based on the object of reducing the problems described above and, In particular, of providing 
an improved estimation apparatus for estimating the transmission operation of a transmission channel for digital data 
which is subject to inference, as well as a receiver and a system for transmitting data using such an estimation appa- 
ratus. 
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[0011] This object is achieved by an estimation apparatus having the features of Claim 1, by a receiver having the 
feature of Claim 12, and by a system for transmitting data having the features of Claim 14. 
[0012] In a particularly advantageous manner, the estimation apparatus according to the invention is used to produce 
at least one estimation sequence which includes characteristics not only of the transmission channel but also of the 
s conversion device, and this results in an improvement in the bit error rate of the overall transmission path comprising 
the transmission channel and the receiving device. 

[0013] It is furthermore advantageous to digital/analogue convert the analogue input signal and to assign a value 
set from a metrics table to the digital/analogue converted value since, in this case, random influences occurring on the 
transmission path can be taken into account using the values from the metrics table. 
io [0014] Those effects of the transmission path which are subject to memory can, furthermore, be taken into account 
in a manner which is efficient in terms of computation, if the values of the metrics table in each case comprise the 
logarithm of the conditional likelihoods for the possible channel states, which results from a finite number of most 
recently transmitted symbols. 

[001 5] In one particularly preferred embodiment, the values in the metrics table become the conditional likelihoods, 
is which in each case cover the likelihood of the occurrence of a received value, subject to the condition that the channel 
is in a specific state, and with this conditional likelihood being obtained by counting the received subsequent values 
from the estimation sequence of estimated channel states. 

[0016] In one advantageous embodiment, a state in each case corresponds to a digital interval of the digital/analogue 
transducer, and an element in the metrics table corresponds to the logarithm of the conditional iikelihood. In conse- 
20 quence, multiplication operations of likefihoods, which are advantageous in terms of data processing, can be expressed 
as the addition of their metrics values. 

[0017] Major advantages result if the conditional likelihoods are obtained by counting the events in that a received 
value occurs when the channel is in a state which is obtained from the estimation sequence (a) since, in this case, 
even a simple memory module can store these values as matrix elements or, in processed form, as elements in a 
25 metrics table. Furthermore, if, after completion of a measurement interval, the likelihood is normalized by forming the 
ratio with the total number of received values which have occurred for each channel state, the pure numerical values 
obtained In this way themselves represent the respective conditional likelihoods. 

[001 8] Dynamic adaptation is advantageously made possible by She elements in the metrics table being updated on 
the basis of the received value sequences. 
so [0019] This concept can be advantageously implemented by means of a data estimator to which the associated 
valueset from the metrics table is supplied and which in each case assigns a most likely value of an estimation sequence 
to the supplied value set and supplies an output signal to a discrete channel estimator which signal describes the 
estimated channel state. 

[0020] In this case, the output signals can advantageously be used for continuous calculation of the metrics table 
36 values in the discrete channel estimator, on the basis of histogram values of the conditional likelihoods. 

[0021] If the time correlations of those elements of interference or nonlinear effects which are subject to memory act 
for a longer time than the duration of only one transmitted symbol, the metrics table values may each be associated 
with more than one previous state change, in this case resulting in corresponding multiplication of the matrix elements 
of the numerical values and of the elements in the metrics table. In this context, reference is also made to the article 
4o by A. Kavcic, J. Moura, The Vrterbi Algorithm and Markov Noise Memory, IEEE Trans. Information Theory, Vol. 46, No. 
1, Jan. 2000 pp. 291-300 whose contents are incorporated completely here by reference. 

[0022] The invention will be described in more detail in the following text using the attached drawings and with ref- . 
erence to preferred embodiments. 



45 Brief description of the figures 



Figure 1 shows a typical optical transmission path in which the invention can be used, 

Figure 2 shows an eye diagram of a signal received from an optical transmission path essentially without 

any interference, in particular the transmission path shown in Figure 1, 

Figures 3a and b show the comparison of the electrical implementation of a conventional detector and of the re- 
ceiver structure according to the invention. 

Figure 4 shows the quantization characteristic of a quantizer, in particular of an analogue/digital trans- 

ducer (ADC), 

Figure 5 shows a block diagram of the assemblies of the receiver structure from Figure 3b, 

Figure 6 shows the transition from numerical values of the respective conditional likelihoods for a quan- 

tized or, in particular, digital value for an estimated channel state to metrics values which indicate, 
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in particular, the logarithm of the respective conditional likelihood, 
Figure 7 shows the relationship between the likelihood density function and the likelihood distribution of 

the quantized or digitized signal. 

s Detailed description of preferred embodiments 

[0024] In the simplest case, an optical transmission system comprises three basic units (illustrated by way of example 
in Figure 1), a transmitter 1 , a transmission medium 2 and a receiver 3. A semiconductor laser whose light intensity is 
modulated with an electrical signal is preferably used for the transmitter 1. 
10 (0025] Various types of optical fibre may be used as the transmission medium 2, in particular such as monomode 
fibres, gradient fibres and the like, which allow propagation guided in the form of waveguides for the electromagnetic 
waves. Commercial systems normally use photodiodes for detection of the optical signals, and these supply a current 
which is largely proportional to the optical power. 

[0026] The task of subsequent detectors is, in general, to reconstruct the original digital data stream 4 with as few 
is errors as possible from the electrical signal 5. 

[0027] By way of example, Figure 2 shows the electrical received signals 5, in each case shifted by an integer multiple 
of a bit interval T and superimposed, of a typical optical.transmtssion path (eye diagram). 

[0028] In conventional systems, the electrical signal 5 Is sampled at equidistant intervals, corresponding to one bit 
interval, and is supplied to a simple threshold-value decision maker 6, illustrated by way of example in Figure 3a. A 
20 "1" or a "0" is detected as part of a signal sequence a depending on whether the signal is above or below the decision 
threshold at the sampling time. 

[0029] The invention is based on the new receiver structure described in the following text and illustrated by way of 
example in Figure 3b. 

[0030] While conventional detectors, as illustrated in Figure 3a, essentially comprise only a low-pass filter 7 and a 
25 threshoid-value decision maker 6, togetherwith a device for samples 8 taken at discrete times, the estimation apparatus 
9 according to the invention provides a receiver structure with lime quantization by means of the sampler 8 and value 
quantization by means of a quantizer 1 0, in particular an A/D transducer 10 and a downstream data estimation appa- 
ratus 11. 

[0031] The electrically preprocessed signal 5 received by the photodiode 3 Is supplied to a filter 12 and Is sampled 
30 by the sampler 8 in time with the symbols. The continuous analogue signal downstream from the filler 12 is converted 
into a digital signal by means of the analogue/digital transducer 10. The output of the A/D transducer 10 supplies the 
data estimation apparatus 11 with a discrete value from a set of a finite number of quantized or digitized values, which 
are associated with the occurrence-of a specific event, for example as the event 1: "the sampled value Is between 
threshold 4" and 5". 

35 [0032] The estimation apparatus according to the invention now tries to reconstruct the transmitted bit sequence 4 
as well as possible from the sampled and quantized values. 

[0033] This object has been achieved by Inventing an estimation apparatus which also includes an estimation device 
for channel estimation 13, see for example Figure 5, which, on the basis of Information arriving from the A/D transducer 
10, adaptively adjusts llself to any desired nonlinear channels, non-Gaussian interference and the natural character- 

40 istics of the A/D transducer 10 or other Imperfections, and compensates for or corrects them as far as possible, and, 
in particular in conjunction with the estimated channel, carries out an estimation of the transmitted bit sequence 4 in 
such a manner that the error likelihood in the estimation sequence a is as small as possible. 
[0034] The invention thus covers two of the previously mentioned problem areas: firstly, static nontinearities and 
dynamic nonlinearities which are dependent on the bit sequence from the A/D transducer 10 are taken into account 

45 and thus compensated for. Secondly, both the resultant conditional inference likelihood densities, which are dependent 
on the bit sequence, and the nonlinear channel model for the probabilistic sequence estimator, are estimated including 
any required unique characteristic, which is necessarily known, of the A/D transducer while data transmission is taking 
place. 

[0035] Based on the simple model of a quantizer which, in this case, is also regarded as the most general form of a 
so digitizer or A/D transducer, one important aspect of the present estimation apparatus will be explained in more detail 
in the following text. 

[0036] The quantizer provides a function which assigns elements in a finite set to the continuous, real input values. 
If r denotes the input value and ? denotes the output value of the quantizer, then 

66 

Equation l q ■ re R fi(r) 3 ?e {? m , ? (2 \ .... r""> 
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[0037] Each input value r e R has only one associated output value 



re {r (U 



, r w 7 W } 



The definition range, which in this case covers the set of real numbers, is therefore subdivided into disjunctive intervals, 

and each interval is assigned an applicable function value r "> . Normally, the number of intervals exactly corresponds 

to the power of the value set. The limits between the intervals are denoted as thresholds. 

[0038] Figure 4 shows an example of a quantization function graphically, in the form of a characteristic. 

[0039] The actual estimation of the transmitted data sequence is based on an apparatus which estimates the state 

sequence of a channel modelled as an automatic system on the basis of theoretical likelihood considerations, and 

uses this to determine the most probably transmitted symbol sequence a. 

[0040] To do this, the data estimator 11 requires so-called metrics increments which are shown, byway of example, 
in Figure 7 as elements "1 " in a metrics table 14 and which are related to the conditional likelihood of a specific analogue/ 
digital converted value and of a channel state, and in consequence are related to the interference likelihood density. 
[0041] One aspect of the present invention is based on obtaining the required metric increments "1 " on the basis of 
channel-state-dependent observation of the information supplied from the quantizer 10. To do this, the events are 
counted as a function of the states, and the overall required statistical characteristics of the channel can be derived 
from this. The metrics increments are stored in a table 14, and the number of events is stored in a store 15, which is 
likewise shown by way of example in Figure 7. Table 14 for the metrics increments "1" is connected directly downstream 
from the A/D transducer 10, as can also be seen in Figure 5. 

[0042] The channel model used In US Patent 5,131,011 generates the possible channel output signals, as a function 
of the state, in the absence of noise, and the channel estimation method also adaptively determines the noise-free 
potential output signals. However, this method cannot be used for optimum consequence estimation for non-Gaussian, 
state-dependent noise. For this reason, the branch metrics Table 14 for the channel model of the invention contains 
the complete stochastic characteristics of the channel from the point of view of the assessment apparatus, which 
characteristics are determined adaptively. 

Theory of an optimum estimator (MLSE) for quantized signals 

[0043] The following text describes the theoretical principles of optimum estimation taking account of an A/D trans- 
ducer and this estimation forms the basis for the inventive solution to the problem of an optimum estimator. 
[0044] The quantizer or A/D transducer provides a function which assigns elements in a finite set to the continuous, 
real input values. If rdenotes the input value of the quantizer and r the output value, then reference is made to equation 
1 , in which case this functional definition includes the requirement for uniqueness, that is to say each input value r s 
R has only a single associated output value 



The definition range {in this case the set of real numbers) is thus subdivided Into disjunctive intervals, and each interval 
is assigned an associated valid function value ? <*> . Normally, the number of intervals just corresponds to the power 
P of the value set. The limits between the intervals [s,.1, sj with which the values r tl! are associated, are referred to 



[0045] If r denotes the vector of the quantized samples, then, on the assumption that successive samples are 
statistically independent, the expression for the maximum likelihood sequence estimator Is: 



r e <.r (1) , 



r< 2 >,.„, r 
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a = arg{nwx(/>{rj*})} 



*{r(n' ( H! 



[0046] Since the random variables can assume only a finite number of discrete values, the conditional likelihood 
density function consists of weighted delta pulses. The a priori likelihood P( r s ]b) of the occurrence of a specific output 
value from the quantizer is therefore used instead of the density function. If it is also assumed that the value r„ sampled 
at the discrete time depends only on the L previously transmitted symbols (finite memory), then: 



Equation 3 

a - arg< 



[0047] Since both the variable f „ and Z Jt can assume only a finite number of different values, the individual branch 
metrics or elements of the metrics Table 14 can store 1( r ^z^) in a table, which is highly conducive to the hardware 
implementation, since this Table 14 can be in the form of a store. 

[0048] A priori likelihood P( ? u|Z M ) can easily be determined exactly from the a priori density function p( ? ? |Z (1 ) of the 
signal with continuous values. In this case: 

Equation 4 



[0049] This means that the a priori likelihood for a specific event corresponds exactly to the area of the a priori 
likelihood density function between the two associated thresholds, as is shown in Figure 7, which explains the rela- 
tionship between the likelihood density function and the likelihood distribution of the quantized signal. 
[0050] In consequence, the variables P( ? B |Z tl ) required for the optimum estimator can, in a practical implementation, 
be obtained by measuring the relative frequency, which corresponds to the number of events related to the total number 
of measurements, with which the samples r w lie in an interval S H <r fl <Sj for a specific channel state Z< k >. However, the 
interval in which a sample lies is precisely the magnitude which a quantizer provides at its output. The digital quantizer 
output can therefore be used as the address of a memory cell, which is incremented by one for the relevant sample 
and the relevant state, and thus contains the number of events. After completion of a measurement cycle, the frequency 
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distribution for each state Z< k > can be deduced from the total number of samples and the memory contents. Since 
likelihoods in logarithmic form are generally used for calculation purposes, the required branch metrics are obtained 
by taking logarithms of the frequency values. 

[0051] The metrics Table 14, which may also be referred to as a branch metrics table, thus contains those metric 
5 increments which can be used directly for a Viterbi algorithm (VA) in the data estimation apparatus 11 . The application 
of the Viterbi algorithm to the values in the metrics Table 14 Is well known to persons skilled in the art, and does not 
require any further explanation at this point. 

[0052] Theoretically, it is also possible to verily that these metrics increments or elements in the metrics table are 
expediently the logarithms of the conditional likelihoods of the respective transitions 1(/,| -- logP(',| zj. This 
10 Table 14 may be implemented by one or more random access memories, in which case the quantized signal then 
represents the respective address. 

Channel estimator 

15 [0053] The channel estimator 13 has the task of estimating the required conditional likelihoods. The channel estimator 
comprises a delay element for the quantized signal r „, a shift register for determining the channel state, and a table 
for counting the events, and the latter is represented as Table 15 in Figure 6. This Table 15 may be implemented by 
one or more random access memories, in which case the quantized and delayed signal and the associated estimated 
state represent the address. The respectively addressed memory cell is incremented by one per symbol clock cycle, 

so so that event counting takes place and a histogram is A thus formed for each state. A 

[0054] if N denotes A the number of events per state z« and it is assumed that H(zO), ? <» )=H(i,j) is the content of the 
memory cell of state zO) and quantized value r l]) , then: 

2fi Equation 5 

z,) = Urn(l/N • H(i,j)), , as N tends to infinity. 

[0055] This means that, provided N is sufficiently large, the event counting carried out results virtually directly in the 
30 conditional likelihood, which is found once again in the expression in the branch metrics For the optimum sequence 
estimator. 

[0056] In order to prevent memory overflows, the memory may be normalized at regular intervals (algorithm, page 5). 
Branch metrics calculation (BMC) 

35 

[0057] The branch metrics are calculated in the block BMC from the histogram H(l,j). To do this, logarithms are taken 
of the individual memory contents H(i,j) in the event store 1 5: 

w Equation 6 

BM(ij) = -ln<H(j,j)) 

[0058] The elements of the metrics Table 14 obtained in this way are then, as is shown by way of example in Figure 
-*s 5, transferred by the discrete channel estimator 13 to the metrics Table 14. A 

[0059] However, it should be remembered that the number of events per state z('> and quantized value may be 
very small or zero. For this reason, adaptation algorithms which are known to persons skilled in the art are used to 
prevent poles from occurring while taking logarithms. 

[0060] it Is furthermore within the scope of the invention to transmit pilot data or test data first of all before signal 
so transmission starts, which data are distributed statistically or randomly over the various channel states by means of a 
scrambler or other suitable devices, and to offer the capability to set the estimation apparatus optimally during a test 
period for it. Such pilot data or test data can, furthermore, be transmitted regularly or as a function of any interference 
that occurs. 



Claims 

1. Estimation apparatus in particular having a channel estimator (13) for estimating characteristics of a transmission 
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channel, in which case the transmission channel is subject to interference, and In which the analogue data (5) 
received from the transmission channel is assigned an estimator sequence (a) of digital data values in particular 
by means of a data estimator (11), 
which is characterized 

5 in that the estimator apparatus produces at least one estimation sequence (a) whose estimation includes both 

characteristics of the transmission channel and characteristics of the receiver. 

2. estimation apparatus according to Claim 1, in which the analogue input signal (5) is preferably filtered and is 
analogue/digital converted as a preferably filtered signal r, and the analogue/digital converted value ? is assigned 

10 a value (ly) from a metrics table (14). 

3. Estimation apparatus according to Claim 2, in which the values in an event table {15) each comprise the number 
of occurrences of a specific analogue/digital converted value in a state (zj of the transmission channel. 

15 4. Estimation apparatus according to Claim 3, in which the state per analogue converted value is obtained from the 
estimation sequence a of the data estimation apparatus (11). 

5. Estimation apparatus according to Claim 4, in which the conditional likelihood P(f ,\[zj of the occurrence of the 
analogue/digital converted value ?» for a state z JC of the transmission channel is obtained from the event table (15). 

so 

6. Estimation apparatus according to Claim 5, in which one element in the metrics table (14) corresponds to the 
logarithm of the conditional likelihood P( 2 f i\zJ of the occurrence of a specific analogue/digital converted value f„ 
in state z^ of the transmission channel. 

25 7. Estimation apparatus according to one of Claims 2 to 5, in which the conditional likelihood Pffplz,,) of occurrence 
of a specific analogue/digital converted value ?, in a state of the transmission channel after completion of a meas- 
urement interval is normalized by forming the ratio with the total number of events N which have occurred for this 
state. 

30 8. Estimation apparatus according to one of Claims 2 to 7, in which the elements of the metrics table (14) are updated 
on the basis of the estimation sequences (a) obtained and the analogue/digital converted values. 

9. Estimation apparatus according to one of Claims 2 to 8, In which tho data estimator (11) Is supplied with the 
associated value set Ijj from the metrics table (14), and which in each case assigns a most likely value of an 

35 estimation sequence (a) to the supplied value set fy and supplies an output signal to the discrete channel estimator 

(13) , from which the channel estimator (13) determines the state of the transmission channel z, r 

10. Estimation apparatus according to Claim 9, in which histogram values are obtained from the elements of the event 
table (15) and are used to carry out continuous calculation of the values of the metrics table (14). 

40 

11. Estimation apparatus according to one of the preceding claims 2 to 10, in which the values in the metrics table 

(14) are each associated with more than one digital/analogue converted value. 

12. Receiver for receiving data transmitted by means of a transmission channel, characterized by an estimation 
45 apparatus according to one of the preceding claims. 

13. Receiver according to Claim 12, further characterized by a filter 12 and a sampler 8. 

14. System for transmitting data characterized by a receiver according to Claim 12 and an estimation apparatus 
eo according to one of the preceding claims 1 to 11. 
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FIG. 3a 



8 

/ 
Sample 



r(t)- 



Low-Pass 
Filter 



0 


















FIG. 4 




_r ^ S3 54 s 5 7q % 
Quantization Characteristic 



10 



EP 1 139 619 A1 



J2 c« 
« E 



+3 ta 

CD Q. 

gZ5 



O 

£ 
2 



< g 

II 



11 



EP1 139 619 A1 



FIG. 6 
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FIG. 7 
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Description 

[0001] The present Invention relates to a method and a receiver according to the preamble parts of claims 1 and 15, 
respectively, in particular the invention relates to the detection of an unreliable detection event (UDE) providing a basis 

5 for error rate estimation. 

[0002] In digital receivers it is often necessary to continuously estim ate receiver performance In terms of the "quality" 
of the detected digital output signal. The on-line estimation of receiver perfo rmance may be used to improve transmis- 
sion quality, for instance by initiating the use of redundant transmission resources providing better transmission quality 
{e.g. channel switching in wireless applications). 

io [0003] Error rate and, more specifically, Bit Error Rate (8ER), is universally accepted as a service oriented 'quality" 
measure of performance of a digital transmission system. Sometimes, an errored-block rate (EBR) is used Instead of 
BER, because the EBR is easier to measure. However, for practical purposes and for not-too-long block lengths, the 
difference between BER and EBR often becomes irrelevant, especially at high signal-to noise ratios (SNR), i.e. at very 
low BER, when most errors occur In isolation. 

15 [0004] Determination of true error rate is clearly impossible without knowing the transmitted bit stream, so one must 
resort to error rate estimation. 

[0005] Unfortunately, In order to minimize error rate, measurements of physical signal parameters that are directly 
related to the controllable receiver parameters, like received signal strength, power, etc., are often Insufficient for con- 
trol, because the dependence of the error rate on these parameters Is unknown, especially for heavily distorted signals 
so which are processed by non-linear processing functions, as is the case e.g. in a maximum-likelihood sequence detector 
(MLSD). 

[0006] Measuring physical signal parameters, such as received signal strength provides only a very crude measure 
of overall system performance In particular the error rate. Hence such measured signal parameters are hardly suitable 
to optimize receiver parameters for minimum error rate. Oniy in some simple applications, such "local" performance 
^ optimization is possible, e.g. when adapting sampling phase towards maximized eye opening in a transmission system 
with hard-decision detection and with a channel without significant Inter-symbol-interference (ISl). But typically such 
optimizations do not work under harder conditions, e.g. when severe IS) occurs or when a more complex detector is. 
used. 

[0007] In a digital communication system symbols are transmitted and typically a number of 2 n symbols are used. 

30 in the binary case (n-1), there are two different symbols, designated logic 0 (zero) and logic 1 (one). 

[0008] In a well-known approach, bit error rate (BER) estimation Is accomplished with the help of a Forward Error 
Correction (FEC) decoder, e.g. by counting the number of corrected bit positions. However, this FEC based estimation 
. does not work well when FEC decoder errors are dominant i.e. at very low SNR. Another drawback Is that the method 
requires a long feedback path from a downstream FEC decoder, which may be designated receiver back-end, to the 

35 receiver front-end, typically comprising demodulator and detector. Yet In some applications no FEC decoder may be 
present, or a framing-agnostic and FEC-agnostlc method of estimation may be used Instead. The latter may be the 
case when the FEC decoder Is administratively separated from the front-end I. e. when the front end is a stand-alone 
product, and hence the receiver back-end may not provide the necessary BER estimation feedback. E.g. the FEC 
decoder may be located In a host system's digital receiver back-end of a separately marketable dense wave division 

40 multiplexer (DWDM) transponder module that comprises a op!ical-to-electrlcal (O/E) and an analog-to-digital (A/D) 
conversion and a symbol detector, but no FEC decoder. 

[0009] In another conventional embodiment, the detected or decoded bit stream Is compared with a known period- 
ically repeated transmission sequence, I.e. a training-based estimation Is carried out. This method has the disad- 
vantage of training sequence overhead, which causes a rate increase, which in turn Is extremely undeslred in case of 
45 bandwidth-limited optical transmission systems, such as DWDM systems, Moreover the method requires knowledge 
of data formats. 

[0010] Bolh FEC based methods and training-sequence-based methods require a framing process and do not work 
on un-synchronized bit slreams. 

[001 1] In a metrics-based estimation approach, decision Information In a Viterbi decoder or detector, which will be 

so explained below, is used for providing a basis of a BER estimation. 

[0012] US 5,944,844, "Method for determining connection quality in a receiver utilizing a Viterbi Decoder*. This in- 
vention uses accumulated palh metric differences, I, d. the sum of decision variables and the minimum of decision 
variables along the decoded path In the VA (Viterbi algorithm) to estimate BER. A drawback of this method is an 
Increase of complexity since decision variables need to be stored until the trace-back step is completed and since 

55 accumulation of decision variable Introduces complexity due to a floating-point format. Apart from the disclosure to 
compare estimates with predetermined threshold values the US 5,944,844 fails to expose a practical mapping from 
the measured observable to a BER estimate. 

[0013] US 6,141,388, "Received signal quality determination method and systems for convolutlonally encoded com- 
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munlcations channels". Unlike the method of US 5,944,844, rather lhan using decision variables accumulated along 
the decoded path, this invention uses only the final decision metfic in a decoded frame. This final decision metric is 
then mapped to a BER estimate. The method is described only in connection with convoiutlonally encoded data streams 
and frame decoding. Again, there Is no disclosure how to map the measured decision metrics to a BER estimate. 

5 [0014] A. J. Viterbi proposed in "Error Bounds for Convotutional Codes and an Asymptotically Optimum Decoding 
Algorithm" {IEEE Trans. Inf. Theory, 1T-13, pages £80 to 269, April 1 967) a decoding algorithm for convotutional codes, 
which has since become known as the Viterbi algorithm. Later on it was recognized that it was in fact a maximum- 
likelihood decoding algorithm for convolutlonal codes (confer Shu Lin, Daniel J. Costello Jr., "Error Control Coding: 
Fundamentals and Applications" Prentice-Hall, inc., Englewood Cliffs, New Jersey 07632, 1983). 

10 [0016] The Viterbi algorithm may also be used for channel equalization In order to cope with ISI. On a binary 1S1 
channel, at a channel memory of m bits, there are 2 m states corresponding to all possible bit sequences of length m 
and 2 transitions entering and leaving each state, i.e. there are 2 m+1 transitions between successive stages or time 
units in the trellis. 

[001 6] In an initializing step of the Viterbi algorithm beginning at an initial stage a path metric for a single path entering 
is each state at the Initial stage Is computed. Each transition between states corresponds to a symbol. The path and Its 
path metric Is stored for each state. 

[0017] The Viterbi algorithm further comprises a repeating step, In the repeating step the path metric for all the paths 
entering a state at a stage is computed by adding the branch metric entering that state to the metric of the connecting 
survivor at the state of the preceding stage, For each state the path with the largest path metric, called survivor path, 

so [s stored together with its path metric, and all other paths are eliminated. 

[0018] A fog-likelihood function log P(rlv) is called the metric associated with the path v and Is denoted M(rlv). The 
different metrics depend on the properties of the transmission path. They may be obtained from noise measurements, 
The metrics may be assumed to be time Invariant and listed In a look-up table for each transition from one state to 
another for a special application, or it may be obtained from on-line measurements {cf. EP 1 139 619 A1 "Channel 

25 estimation using a probability mapping table"). 

{0019] in a truncation mode Viterbi detector, path memories are truncated after D symbols and a decision is 
enforced on transitions to the states that are leaving the "window" of D+-1 Irellis stages. More specifically, such a detector 
has a path memory that contains D 1 1 trellis stages resulting In D transitions, When computing survivors from time t+1 , 
a decision on the detector output is forced for a specific transition from a state at time instant t-D to a state at time 

so Instant t-D+1, As usual, in a practical Implementation the states on this selected transition could be the trace-back 
states on the best path at time t, or they could be traced back from any state at lime t. 

[0020] A Viterbi detector operating in block-mode (block-mode VA) detects bits in a block i.e. a trellis segment by 
finding the best path between two previously determined states at the block boundaries. For the purpose of this Inven- 
tion, it Is immaterial which detailed method Is used to determine the states at the block boundaries. 
35 [0021] For highest throughput applications like high-speed optical communication systems operating at data rates 
e. g. at 10 Gbps or 40 Gbps or above, It is essential that the error rate estimator can be Implemented with very low 
additional complexity. 

[0022] It Is the object of this Invention to provide a method and a receiver using another observable which does not 
require extensive calculations and/or additional circuitry. 

40 [0023] This object Is achieved by the subject matter of the Independent claims. 

[0024] Preferred embodiments of the invention are the subject matter of the dependent claims. 
[0025] This invention advantageously does no! require a framing process and can operate on non-synchronized bit 
streams, unlike both the FEC-based methods and the training-sequence-based methods. So the present invention 
does not require a FEC code nor a known training sequence for error rate estimation. 

45 [0026] Advantageously the same principle can be applied for symbol detectors, maxlmum-llkellhood-sequence de- 
tectors, truncation mode sequence detectors, block processing sequence detectors, and frame decoders for trellis 
encoded frames. This means that this Invention can be applied to all detectors and decoders where a random decision 
can be Identified. 

[0027] In the binary case I.e. if only two symbols designated logic 0 and logic 1 are sent and detected on the receiver 
so side the unreliable detection event can be determined by checking as to whether the analog voltage Is within a range 
around the threshold which is used for discriminating logic 0 and logic 1 . The tatter step is a kind of an analog-to-digital 
conversion. 

[0028] As a further advantage, the additional condition may be used that a decision must be a relevant random 
decision, i.e. a random decision where a wrong decision necessarily leads to a bit error. This additional condition 
55 ensures that only events are considered which are related to the bit error rate, by definition. 

[0029] The advantage of considering more than one unreliable detection events by different thresholds for the ab- 
solute value of the metric difference between survivor and competitor generates more observabtes for the error rate 
e.g. the unreliable detection events obtained by using a lower threshold may be allowed for In the case of lower bit 
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error rates and vice versa. 

[0030] The condition that two metrics are equal in case that the metrics are represented by digital values can be 
verified by a single operation in most processors after the metrics have been loaded Into appropriate registers. This 
Illustrates that dedicated logic for checking the equality of two registers requires only a small amount of dedicated 
s circuitry. 

[0031] Counting unreliable detection events during a fixed period of time and reading the counter value thereafter 
results in a value which Is proportional to the frequency of the unreliable detection events, which In turn Is related to a 
bit error rate approximately by simple equation (8). 

[0032] The adaptation of the period of lime for counting the unreliable detection events to the frequency of the un- 
10 reliable detection events may ensure quick adaptation of the receiver parameters In orderto minimize the bit error rate 
and to ensure a sufficiently smail uncertainty In the frequency of detection events. 

[0033] If thecountervaluelsconsldered to betoosmalli.e. lowerthanafirsf threshold, thecountervalueispreferrably 
not reset, such that the current measurement period is extended and accumulated values are not discarded.. 
[0034] The setting of demodulator parameters In order to minimize the unreliable detection event frequency and 
is thereby the bft error rate keeps the bit error rate at an optimum even If the channel properties change In time. 

[003S] In the following preferred embodiments of this Invention are described referring to the accompanying drawings. 
In the drawings: 

Fig. 1 shows an optical transmission system, 

20 

Fig. 2 Illustrates the unreliable detection event for binary symbols; and 
Fig. 3 and Fig. 4 show simulation results. 





Abbreviations: 


A/O: analog-to-digital 


ADC: analog-to-digltal converter 




AGC: automatic gain control 


MLSD: maximum-likelihood sequence 


so 


AWGN: additive white Gaussian noise 


detection 




BER: bit error rate 


O/E: optlcal-to-electrical 




DWDM: dense wave division 


PDF: probability density function 




multiplexer 


PMD: polarization mode dispersion 




EBR: errored-block rate 


RDE: reliable detection event 


35 


FEC: forward error correction 


SNR; signal-to-noise ratio 




ISI: Inter-symbol Interference 


UDE: unreliable detection event 






VA: Vfrerbl algorithm 



[0036] An embodiment comprises a low-complexity error rate estimation method based on information that Is readily 
available In a Maximum-Likelihood-Sequence-Detector (MLSD). It Is especially useful for a local error-rate-based con- 
trol loop within a digital receiver front-end (demodulator and detector), 

[0037] More specifically, detection performance of a MLSD receiver, in terms of minimizing post-detection error rate 
is optimized, by providing a low-complexity post-detection error rate estimate to a receiver control facility for unframed 
data streams. 

[0038] To this end "unreliable detection events" (UDE) are detected, counted, optionally processed, and then used 
to estimate error rate. 

[0039] Fig. 1 shows an optical transmission system. It comprises a transmitter 1 , and optical link 4 and a receiver 5. 
A typical transmitter 1 comprises an FEC encoder 2 for encoding input data a, In order to generate data d,. This data 
d. Is forwarded to a modulator 3. The modulator 3 generates an optical sig na! y(t) constituting the output of transmitter 1 . 
[0040] The optica! signal is transmitted via optical link to receiver 5. The optical link comprises optical fibers which 
attenuate the optical signal and in addition constitute a dispersive channel. In order to compensate for the attenuation 
the optical link may comprise optical ampliliers which add noise to the optical signal. In state of the art dense wave 
division multiplexer (DWDM) systems the optical signal suffers from severe signal distortions that are caused by chro- 
matic dispersion orgroup velocity dispersion, polarization mode dispersion (PMD), self-phase-modulation, inter-symbol 
interference {ISI) caused by the receiver itself, and other signal distortions that can beneficially be removed or coun- 
teracted by means of a MLSD. 

[0041] At the receiver side of the optical link the optical signal r{t) is Input into receiver 5. Receiver 5 comprises 
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physical interface 10 which performs an O/E conversion. The analog electrical signal may be directly input into auto- 
matic gain control (AGC) circuit 20. Both, the physical interface 10 and AGC circuit 20 comprise an upper cut-off 
frequency. Both cut-off frequencies must be higher than 1/{2T), T being the symbol period. On the other hand too much 
excess bandwidth in excess over the required minimum picks up more noise from the optical link 4, which degrades 

5 the receiver performance by increasing the bit error rate. In typical receiver designs an excess bandwidth of 50% to 
100% is therefore provided {S. U. H. Qureshl, "Adaptive equalization", Proc. IEEE, Vol. 73, 1 985, pp. 1349-1387). 
[0042] In practice due to the design complexity of low pass filtering in the GHz range, often only implicit low pass 
filtering is performed In the receiver, by virtue of given component bandwidth limitations. However, In order to obtain 
a more precise control of the bandwidth, another embodiment may comprise an adjustable low pass filter 18 between 

10 physical Interface 10 and AGC circuit 20. The low pass filter 1 5, especially its cut-off frequency, Is controlled by output 
86 of control unit 70 in order to minimize the bit error rate. 

[0043] The AGC circuit 20 may amplify the analog electrical signal output either by physical interface 1 0 or low pass 
filter 15 to a constant level In terms of average rectified voltage or root-mean-square voltage. In another embodiment 
the amplification of AGC circuit 20 may be controlled by control unit 70, output 85 based on digitized values rj (cf. US 
is 3,931,584). 

[0044] The sampling phase of A/D converter (ADC) 40 may be controlled by a sample-and-hold circuit 30 illustrated 
by a switch In Fig. 1 . The sample-and-hold circuit 30 is controlled by output 84 of control unit 70. The ADC 40 performs 
a three-bit conversion and outputs digital data rj. Various parameters of the ADC 40 may be controlled by control unit 
70. Therefor outputs 81 to 83 of control unit 70 are provided, Output 81 may control the offset of ADC 40, output 82 
so may provide a reference voltage which controls the amplification of the A/D conversion I. e. the spacing of the thresholds 
and output 83 may conlrol the non-linearity of the ADC 40. 

[0045] In other embodiments an N-bit conversion may be performed by the ADC 40, wherein N is an integral positive 
number. 

[0046] The digitization of the analog signal output by AGC circuit 20 may be performed at the symbol rate of 1/T. In 
25 another embodiment n-fold over-sampling may be performed i.e. the analog voltage output by AGC circuit 20 Is sampled 
at a rate of 1/(nT). The sampling rate is controlled by control unit 70, which controls the sample and hold circuit 30 and 
the digitization of ADC 40. In high-speed applications n Is often chosen to be 2 since higher values of n are not feasible 
due to component restrictions. In the case of n=2 the rising and falling clock edges can be used to trigger sampling 
based on a recovered symbol clock. The symbol tj output by ADC 40 Is Input into MLSD 50 which may be Implemented 
so by a Vlterbl decoder. More specifically the Vlterbl decoder may be of a truncation mode type or a block mode type. As 
explained abovo a Vitcrbi decoder determines the most likely path, designated survivor, on the basis of metrics. The 
metric of the most likely path Is output at output of MLSD 50. The metric of the best competitor path is output at 
output m to of MLSD 50. 

[0047] in the case of a Truncation Mode Vlterbl Detector a UDE is declared if the survivor decision at time t-D was 
3S a random decision. The only addition to the standard Viterbt detector is an additional memory containing the Information 
about path metric equality, one bit per stage and state. This memory Is required due to the delay of D steps between 
detecting path memory equality and using that decision. 

[0048] In another embodiment a UDE can also be defined and detected if a VA Is run to determine a most likely slate 
at a particular stage of the trellis, starting from identical zero path metrics. A UDE may only be declared If this best 
40 state selection is subsequently used for detection, i.e. when it Is relevant. 

[0049] In some trellises estimation of bit error rate can be improved if a UDE is declared not for arbitrary path metric 
equalities, but only for path metric equalities between the best path that detects a one on a specific transition and the 
best path that detects a zero on that specific transition. 

[0050] MLSD 50 outputs the detected symbols of the most likely path X, to FEC decoder 60. The FEC decoder 60 
« finally outputs data Uj. In one embodiment the FEC decoder 60 is designed to work with a relatively high Input BERof 
e.g. 10" 4 and still provide a very low output BER of e.g. 10-1 5 , This situation is common In modern optical high-speed 
systems such as DWDM systems employing so-called advanced FEC. 

[0051] Unreliable detection eventsmay be declared If themeirics of the survivor and competitor path are equal which 
Is detected in comparator 51 . A second kind of UDE may be declared If the absolute value of the difference between 

so the metrics of the survivor and competitor paths are lower than a threshold t. This condition Is checked in comparators 
52 and 53 for two thresholds t t and t 2 . There are embodiments of this invention that comprise only comparator 51 or 
comparator 52. Other embodiments comprise both comparators 51 and 52. In further embodiments more than two 
comparators 52 and 53 maybe provided with or without comparator 51, wherein each of the comparators 52 and 53 
checks for a different threshold t. 

55 [0052] MLSD 50 provides a decision clock on line 57 to counters 54, 55 and 56. The decision clock comprises a 
pulse each time a decision is to be evaluated. In the case of truncation mode Viterbi detection, the MLSD generates 
a decision clock pulse each time a decision is enforced after D+1 trellis stages. In the case of block mode Viterbi 
detection, a decision clock pulse is generated each time a block of data Is detected. Each of the counters Is enabled 
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by the decision clock on line 57 and is incremented if the respective comparator has detected an UDE. The counters 
are read from time to time by control unit 70 via lines 87 to 89 and reset if decided by control 70. The counter values 
are proportional to the frequencies of unreliable detection events that are evaluated In control unit 70. The control unit 
70 may read and reset the counters after a fixed counting period. This period may be different for each counter. 
s [0053] In another embodiment the counting period of each counter may be adapted to the frequency of UDE. To this 
end two thresholds may be provided for each counter. 

[00S4] If the counter value is lower than the first threshold, the counting period Is considered to be too short. Con- 
sequently the counting period is increased and the counter may be reset. Alternatively the counter value is not reset 
and counting Is continued for the difference between the old and the new counting period before the counter value is 

10 read out again and compared to the two thresholds, In high frequency applications, the operations of reading out the 
counter value, comparing It to two thresholds, calculating a new counting period and the difference between new and 
o!d counting periods can not be completed within one detection clock cycle. Consequently one ore more UDEs may 
be detected and counted meanwhile. In order to avoid Incorrect measurements the counter may be disabled until the 
new counting period is determined. 

is [0055] if the read out counter value is greater or equal to the first threshold and smaller or equal to the second 
threshold, the counter value is considered to be In the right range and the counting period remains unchanged. After 
reading and processing the counter value the respective cou nter is reset. 

[0056] It the counter value Is greater than the second threshold (he counter vaiue Is considered to be too great and 
consequently the counting period Is decreased. After reading and processing the counter value the respective counter 
so is reset, 

[0057] In order to avoid too many changes of the counting periods and the computational effort associated herewith, 
the counting periods may be kept constant for a predetermined number of counting periods. The latter additional rule 
may only be applied afterthe first time after power up or reset the con ntlng period is within the first and second thresholds 
i. e. has reached the target range. 
S5 [0058] The control unit 70 may average or smooth read out counter values or the frequencies of UDEs obtained 
therefrom. The control unit can obtain BER estimates from the UDE frequencies or counter values by equation (8), 
which will be explained below. 

[0059] The control unit 70 may perform a gradient search over the estimated error rate In order to set one or more 
receiver parameters In order to arrive at an error rate as low as possible. More specifically, in the absence of 'locally" 

30 measurable "physical" quality indicators, steepest-descent minimization of estimated error rate can be used to optimize 
such parameter settings, at least when the error rate has a convex functional dependency on the parameter(s) in 
question, I.e. where a unique globally minimal error rate exists over the relevant parameter space. Error rate optimi- 
zation is especially beneficial if channel parameters e.g. properties of optical fiber 4 or physical Interface 10, often 
comprising a photo diode, change gradually. In other embodiments the control unit may perform the gradient search 

35 on one or more UDE frequencies or counter values. Moreover the control unit 70 may set the gain of the automatic 
gain control circuit 20 via output 85. A further receiver parameter is the cut-off frequency of low-pass filter 15 which 
can be set via output 86 of control unit 70. Further the control unit 70 can control the sampling times by sample-and- 
hold circuit 30 via line 84. Then the control unit influences the thresholds of the ADC 40 via outputs 81 to 83. 
[0060] Skilled persons will understand that not all circuit parts may be present In receivers according to this Invention . 

40 Controlling the sampling phase In sample and hold circuit 30 Is especially beneficial In receivers that do not perform 
over-sampling. 

[0061] If a high order low-pass filter Is used for low-pass filter 15 the high order low-pass filter generates additional 
Inter-symbol Interference, which can be handled by the MLSD 50, There Is a performance trade-off between noise 
reduction and Increased ISI, which varies with channel parameters. The control unit 70 can adjust the receiver param- 
os eters to gradually changing channel parameters by a gradient search over e.g. estimated error rate as explained above. 
[0062] In the following the concept of unreliable detection events will be explained in more detail. 
[0063] A relevant decision Is a decision made in the detector, which is part of a receiver, such that a wrong decision 
necessarily leads to one or more errors in the detected output symbol stream. For Instance, in a Viterbl Detector, 
selecting the wrong incoming branch as a survivor at some state in the trellis Is not a relevant decision If the selected 
so branch eventually Is not on the survivor. Since the MLSD 50 outputs a decision clock at line 57 Indicating a final decision 
on the path all decisions taken by comparators 51 , 52 and 53 are relevant decisions. 

[0064] A decision is called a random decision if it is made without sufficient information which alternative is more 
likely to be correct than the otherfs). A strict random decision Is a random decision under total indifference amongst 
alternativefs), e.g. when a decision is made in a Viterbl detector between two alternative paths with exactly Identical 
55 metrics. This means that decisions detected by comparator 51 are strict random decisions. 

[0065] A wide-sense random decision is a random decision under approximate indifference amongst alternative^), 
e.g, when a decision is made in a Viterbl detectorbetween two alternative paths with an absolute path metric difference 
that is smaller than some small, predetermined threshold value. Consequently comparators 52 and 53 detect wide- 
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sense random decisions. 

[0066] In this patent, an Unreliable Detection Event (UDE) can be defined as an event that fulfills one or more of the 
following conditions: 

5 • The event occurs in a detector that is part of a receiver. 

» The detector makes a random decision. 

• The random decision is a relevant decision, i.e. it affects the detected output symbol stream. 
10 The events detected in comparators 51 , 52 and 53 fulfill all of the three conditions above. 

Instead of an exact theoretical derivation of the relationship between the UDE frequency and 8ER or EBR two 
simplified cases are considered, 

The first scenario is "artificial 0 and only motivates, why observing the frequency of unreliable detection event 
can give information about the true error rate, In a Idealized setting. The results of this scenario apply to both, 
is Symbol or Sequence Detection. 

A hypothetical binary "genie" receiver R with the following properties Is considered: 

[0067] Receiver R decides either compietely reliable I, e. error-free or completely unreliable i. e. at random, for a 
so symbol or for a sequence e. g. a block, and R also knows perfectly when an unreliable detection event has occurred 
because it provides perfect reliability Information. 
[0068] It is assumed that 

• the probability of the event UDE to detect unreliably is P ud 

SB 

• The probability of the reliable detection event (RDE) to detect reliably Is J-P Md 
[0069] In this case 

SO » the conditional probability P(elRDE) of the event e to cause an error, In case of rcllablo detection, is 0.0 (perfect 
decision) 

• the conditional probability P(elUDE) of the event e to cause an error, in case of unreliable detection, Is 0.5 (random 
decision) 

35 

[0070] The total probability of error P(e) at the output of R is P(e) = 0.0 P ud + 0.5 P ud which provides a linear relation 
between the probability P ud of an unreliable detection event and the error probability being equivalent to the error rate. 
[0071] Hence, an observer counting the occurrence of unreliable detections, as declared by R, over a sufficiently 
period of time, can estimate P ud from the relative frequency of unreliable detection events. By virtue of the relation 

•to BER = 0.5 P ud , he can then estimale the error rate. 

[0072] The second scenario is for a binary symbol detector on Additive White Gaussian Noise (AWGN). This case 
Is much easier to model mathematically since it is simpler than a MLSD as addressed in this Invention. 
[0073] The method described there could be generalized as another BER estimation method for symbol detectors. 
Actually, it is a special case of the general approach to derive an estimate of the error rate from a measurement of the 

4s amplitude distribution function. Knowing the amplitude distribution function empirically, model distributions for the con- 
ditional distributions for each possible symbol could be fitted, and the BER could then be derived analytically from 
these distributions. Similar techniques are used in simulation software such as VPItransmissjonMaker™, albeit based 
on full knowledge of the transmitted bit sequence. 

[0074] Assume a binary random data stream transmitted using binary amplitude modulation over a memory-fess 
so channel with AWGN and a bit detector. The a-prlori probabilities of sending a one, P(1), and of sending a zero, P(0), 
are assumed to be idenlical. Further noise Is assum ed to be sig nal-independent, I.e. identical additive noise distribution 
and hence Identical variance a a for received one and zero, More specifically, let u, (104) and u 0 (1 03) be the amplitudes 
of a received one or zero, respectively, In the absence of noise (cf. Fig. 2). Under Ihe above assumptions the conditional 
probability density function (PDF) of observing amplitude u, conditioned on the event of a sent one or zero, is 
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P („|,) ' exf^l (1) 

and 

p (u\0^'J-^^\ (2) 
respectively (cf. 108, 107, respectively). 

[0075] The optimum decision threshold Is at the amplitude where the two conditional PDFs assume equal values, . 
which is at the threshold 0=(u 1 +u o )/2 (cf. 105). The graphs ofp(iiM) and p(ulO) as well as 0 are plotted in Fig. 2. 
[0076] With this decision threshold, the observed probability of bit error P b (e), alias BEfi, Is known as 

where d is the distance between u s and u 0 , d=lu 1 -u 0 l. With the well-known approximation 



35 (for x -> «), we can write this as 



[0077] On the other hand, unreliable detection occurs when we observe a received amplitude u with approximately 
equal distance from Uj and u 0 Is observed, I.e. when 0=u, In this case lu-u.|Wu-u 0 !=d/2 and equation (5) Is obtained 
45 from equations (1) and (2). 



[0078] In this case, deciding for either one or zero is a strict random decision, i.e. the decision is incorrect with a 
probability of 0.5. 

[0079] The probability mass P^ of an unreliable detection event P(0 - Au 12 <, u 0 + Au / 2) Is then approximately 
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m Ah • (P{\)-p(u = 0 1 \)+P(0)- p(u = 0 1 0))= A« • 




(6) 



[0080] Taking logarithms and choosing Au=y-d as a multiple of d, equations (6) and (7) are obtained: 
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in ^)4(f) 2+, {!) +,n(Y) - |n( ^ 



(6) 



ln( 



i(Pb(e))-ln(P ud )+ln(2)-ln( Y )~2lr 




(7) 



so 



[0081] For large d/a, the log-probability ln(P ud ) Is dominated by the term quadratic in d/a. and ln(P ud ) equals ln(P b 
(e)) apart from additive constants and apart from the term logarithmic In d/a which becomes negligible for large d/o, I. 
e, for low error rates. 

[0082] The relationship between logarithmic bit error probability and logarithmic probability of an unreliable detection 
25 event Is monotonlc up to some high BER value, for d/o > 1 . Hence it Is Invertible In this range. In general, P ud can thus 
be used to predict P b (e) approximately provided that y is known, or it can be mapped via a lookup-table. 
[0083] For instance, in the context of sampling phase adjustment for a given channel, noise variance a 2 may be 
assumed to bo constant, but the distance d between mean values of one and zero will vary, and consequently BER 
will vary as well. When we assume that an AGG circuit Is used, d will remain about constant but the noise o e will be 
SO amplified such that again tho "signal-to-nolse ratio" <Pla 2 will vary. 

[0004] For y=1/8 It can be verified, that P b (e) in this case provides an upper bound of true 
[0085] BER in tho range of BER from 0 to 3- 1 0" 3 . This choice of ^corresponds roughly to the case of 3-bit-qu animation 
with automatic gain control. The A/D conversion thresholds for a 3-bit conversion are designated 101 In Fig. 2. The 
digital values are designated 102. 
35 [0086] For practical purposes equation (7) may be simplified to: 



4° [0087] A may be considered to be constant for low error rates. As can be learned from Fig. 4, Equation (8) may be 
used not only for symbol detectors but also for sequence detectors. In the simulation the results of which are shown 
In Fig. 4, equation (8) gives a reasonable estimation for the BER in the whole range of the BER from 3'1Q' 4 to 3-10 1 , 
It Is to be expected that the value of A would change for a different parameters e.g. SNR. Moreover there might be 
circumstances in which equation (8) can not be used. 

45 [0088] Fig. 3 and 4 show simulation results of a transmission of about one million bits. The sampling phase was 
changed in an MLSD receiver. Using about 10 6 bits implies that, in this simulation setup, a BER below 10' 5 cannot 
estimated meaningfully because the observation Interval Is loo short. 

[0089] These simulations have been performed for optically amplified links i.e. optical signal-to-nolse ratio limited 
transmission for the case of severe chromatic dispersion, I s1 order PMD and significant IS1. 
so [0090] While the present invention Is described with reference to the embodiments as Illustrated In the following 
detailed description as well as in the drawings, It should be understood that the following detailed description as well 
as the drawings are not intended to limit the present invention to tho particular Illustrative embodiments disclosed, but 
rather the described illustrative embodiments merely exemplify the various aspects of the present invention, the scope 
of which Is defined by the appended claims. 
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Claims 

1. A method for a receiver comprising: 

s receiving a seq uence of symbols; 

characterized in that 

determining (51 , 52, 53) an unreliable defection event at whtch the absolute value of a difference between a first 
metric (m ls ) for a first event minus asecond metric m lc for a second event Is lower than a predetermined threshold, 

10 

2. The method of claim 1, characterized in that said first metric monotonously depends on the likelihood that a 
received symbol (u, 1 06) said sequence of symbols is a first symbol (u 0 , 1 03) and said second metric monotonously 
depends on the likelihood that said received symbol (u, 106) Is a second symbol (u f , 104). 

3. The method of claim 2, characterized in that said first symbol is a logic "0" and said second symbol Is a logic "1°, 
said unreliable detection event Is determined If said received symbol is within a first range {102); said first range 
comprising a threshold (Q, 105) used for discriminating as to whether said received symbol is said first or said 
second symbol. 

2" 4. The method of claim 1 , characterized In that said sequence of symbols being decoded by a maximum-likelihood- 
sequence detector (50); said first metric (m, s ) being a path metric associated with the most likely path designated 
survivor and the second metric (m fc ) being a metric associated with the second most likely path, designated com- 
petitor, 

25 5. The method of one of the claims above, characterized in that an unreliable detection event Is only determined If 
a decision Is based on said first and second metric and if said decision necessarily leads to at least one bit error 
if said decision Is wrong, 

6. The melhod of claim 4, wherein said maxlmum-likellhood-sequence detector is a truncation-mode Viterbi detector 
30 which enforces a decision after a predetermined number of symbols. 

7. The method of claim 4, wherein said maximum-likelihood-sequence detector is a block mode Viterbi detector, 
wherein the survivor Is the finally decided path in a block and tho second metric being the metric associated with 
the second most likely path In said block. 

35 

8. The method of one of claims 1 to 7, wherein a first unreliable detection event is defined by a first predetermined 
threshold (51) for the absolute metric difference and a second unreliable detection event is defined by a second 
predetermined threshold (52, 53) for the absolute metric difference. 

■to g. The method of one of claims 1 to 8, wherein said first and second metrics are represented by digital values and 
said first and second metrics are equal (51 ). 

10. The method of one of the preceding claims further comprising: 

« counting the unreliable detection events by a counter (54, 55, 56); 

reading the counter value after a first period of time. 

11. The method ol claim 10, further comprising: 

so 

increasing said first period of time, If said counter value at the end of said first period of time is lower than a 
first threshold; 

reselling said counter and maintaining said first period of time, if said counter value at the end of said first 
ss period of time is greater than said first threshold and smaller than a second threshold; and 

resetting said counter value and decreasing said first period of time, ii said counter value at the end of said 
first period of time Is greater than said second threshold. 
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12. The method of claim 11, further comprising: 

setting a demodulator parameter (20, 30, 40} depending on the value read from said counter on the basis of 
a gradient search. 

5 

13. The method of claim 12, characterized In that said demodulator parameter Is an amplification (20) within the 
receiver before an analog-to-dlgital conversion (40). 

14. The method of claims 12 or 13, characterized In that said demodulator parameter Is the sampling phase (30), a 
io quantizer bias (81 ), a non-linearity parameter (83) , an u pp er cut-off f requency of a low-pass fitter (15, 86) a residua! 

chromatic dispersion compensation or a phase association. 

15. A receiver comprising: 

15 a metric generator (50) for generating a first metric (m b ) for a first event and a second metric (rry for a second 

event, 

a unreliability detector (51, 52, 53) for detecting an unreliable detection event when the absolute value of a 
difference between said first metric minus said second metric Is lower than a predetermined threshold. 
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FIG.2 
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Description 

[0001] This Invention pertains to a self-timing method for phase adjustment of a sampling phase according to the 
preamble part of claim 1 and a digitizing, self-timing circuit according to the preamble part of claim 7. More specifically 

s this Invention relates to a method and a circuit for adjusting the phase of a symbol cloGk recovered from an analog 
signal that does not provide a carrier in order to minimize the bit error rate after symbol detection. 
[0002] A method of the preamble part of claim 1 and a circuit of the preamble part of claim 7 are described in WO 
02/30035 A1 titled "Symbol Timing Recovery Method for Low Resolution Multiple Amplitude Signals'. 
[0003] Modern receivers are used to demodulate, detect and decode digitally transmitted signals e. g. in mobile or 

10 fiber optical communications. Due to varying power level of the Input signal, an Automatic Gain Control (AGO) circuit 
Is normally Included before an Analog-to-Dlgltal Converter (ADC) in order to adapt the signal amplitude to" the conver- 
sion range. A digital equalizer (EQ) and an error correcting circuit (ECC) usually follow the ADC in order to arrive at 
an acceptable bit error rate (BER). A clock recovery (CR) circuit extracts frequency and phase and generates a local 
sampling clock at approximately equidistant points in time, with some more or less fixed phase relation to the transmitted 

t« symbol stream. A Sampling Phase Adjustment (SPA) circuit may or may not be present. Its purpose Is to add a phase 
delay to the recovered sampling clock, in order to oplimize receiver performance. In simple receivers for low-distortion 
channels, the SPA is usually missing. 

[0004] Various clock recovery schemes are discussed in 'integrated Fiber-Optic Receivers" by A. Buchwald, K. W. 
Martin, Kluwer Academic Publishers, 1995 (later referredto as B Bucfiwald95°). The US 5,048,060 titled "Digital Signal 
so Receiving Circuit with Means for Controlling a Baud Rate Sampling Phase by a Power of Sampled Signals' describes 
a clock recovery method according to which an analog input signal Is sampled to produce a series of sampled signals 
and adjusts the sampling phase based on the series of sampled signals. To this end the sampled signals are equalized 
and digitally filtered. 

[0005] In practlca! digital commun Icatton systems for high-speed transmission and In mobile communications through 
2f a band-llmlted channel, the frequency response of the channel Is not known with sufficient precision. Moreover, many 
channels have a priori unknown and/or time-variant frequency-response characteristics, like, e.g., wireless channels 
or optical channels. For such channels, it Is not possible to design optimum fixed demodulation filters; they require 
adaptive equalization, instead. 

[0006] Channel distortion results In Intersymbol Interference (IS1), which, when uncompensated, causes high error 
30 rate. This calls for receivers v/lth equalizers that can compensate or reduce the IS! in the received signal, There are 
three well-known types of equalizers: linear equalizer, decision feedback equalizer and maximum-likelihood sequence 
detector also known as Vilerbi detectors (cf. S. U. H. Qureshl, "Adaptive Equalization", Proc, IEEE, vol. 73, 9, pp. 
1 349-1 387, 1 985; later referred to as "Qureshi 85"). 

[0007] In symbol rate equalizer structures, (he equalizer taps are spaced at the reciprocal of the symbol rate. 

35 When the channel characteristics are unknown, the receiver filter Is sometimes matched to the transmitted signal pulse 
and the sampling time Is optimized for this sub-optimal filter. In general, this approach leads to equalizer performance 
that Is very sensitive to the choice of sampling lime (cf. J. G. Proakis, "Digital communications", McGraw-Hill, 2001). 
[0008] In contrast to the symbol rate equalizer, a fractionally spaced equalizer (FSE) Is based on sampling the 
incoming signal at least as fast as ihe Nyqulst rate. in general, a digitally Implemented fractionally spaced equalizer 

40 has tap spacing of mlJn where m and n are integers and n>m, and T Is the symbol period. Usually, a T/2-spaced 
equalizer is used in many applications. Simulations results demonstrating the effectiveness of the FSE over a symbol 
rate equalizer have been given In Qureshi 85; G. Ungcrboeck, "Fractional tap-spacing equalizer and consequences 
for clock recovery In data modems", IEEE Trans. Commun., vol. COM-24, pp. 856-864, 1976 and S. U. H. Qureshi 
and G. D, Jr., Forney, "Performance and properties of aT/2 equalizer", Natl. Telecom, Conf. Record, pp. 11.1.1-11.1 .14, 

45 Los Angeles, 1977. The authors have shown that the FSE does not exhibit sensitivity to choice of timing phase. How- 
ever, for severely distorted channels this conclusion is not valid. Hence, timing recovery means for providing an ap- 
propriate sampling phase have to be Implemented even for a FSE. 

[0009] Under the crucial assumption that the received signals are sampled at proper sampling freq uency and phase, 
many techniques are known that are capable of accurately restoring the received signals into original digital signals. 

bo Therefore, one of the most Important processes in a digital receiver Is to accurately restore symbol timing in terms of 
frequency and phase, i.e. to delermine the poinls in time when Ihe received signal should be sampled. In some older 
systems, a specific carrier corresponding to a sampling freq uency was transmitted, providing frequency and phase for 
sampling, However, most modern systems do not provide a carrier; It Is therefore important to accurately recover a 
symbol timing designated also clock recovery (CR) using only the received signal designated self-timing. Without ex- 

55 pliclt frequency and phase information being available in Ihe receiver, receiver complexity is increased e.g. due to clock 
recovery circuits or due to the need for error resilient facilities, like error correction circuits, 
[0010] Generally, digital methods for sampling phase adjustment can be divided into two groups: symbol-rate sam- 
pling methods and oversampling methods utilizing a clock, which Is faster than symbol rate. 
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[0011] In Symbol-rate Sampling methods, data Is sampled at symbol rate in order to determine sampling phase 
based on optimization of some control parameter. 

[0012] In US 4,494,242 titled "Timing Recovery In a Baud- Rate Sampled-Data System 0 , timing recovery Is based on 
the recognition that the precursor portion ot the response of the channel is relatively Invariant and predictable. By 
selecting a specified threshold point In a small-amplitude region of the precursor portion, a basis Is provided for gen- 
erating accurate bipolar error signals that insure consistent baud-rate sampling of received pulses. A similar approach 
is disclosed in US 4,959,845 titled "Receiver of a System for Transmitting Data Symbols at a Given Baud-Rate". 
[0013] However, the assumption of an invariant and a predictable pulse shape fails in many time-variant channels. 
In US 5,048,060 titled "Digital Signal Receiving Circuit with Means for Controlling a Baud Rate Sampling Phase by a 
Power of Sampled Signals", Inventors calculate an element of the autocorrelation function of a series of baud rate 
sampled signals and use its value to find the correct sampling phase, Actually, they choose the sampling phase that 
either maximizes or minimizes this element of the autocorrelation function. This approach is not applicable in severely 
distorted channels such as, e.g., optical channels with substantial chromatic dispersion and polarization mode disper- 
sion. For example, It can be shown by simulation results that for special channel conditions the autocorrelation function 
values R{0) and RfT), where T Is a bit period, vary very little with sampling phase. Consequently finding extrema is 
actually Impractical. For other channel conditions, the rolovant oxlremum may even occur at a sampling phase with 
non-minimal BER. Also, if an (AGC) circuit uses outputs of an ADC for amplification adjustment, the AGC circuit will 
tend to keep signal power at a constant value and the power measurement cannot be used for timing recovery. 
[0014] An Interesting method is described in WO 00/19655 titled "Timing Recovery for a High Speed Digital Data 
Communication System based on Adaptive Equalizer Impulse Response Characteristics 0 . The timing recovery tech- 
nique is based upon an analysis of the impulse response of the equalizer structure having a feed forward equalizer 
and a decision feedback equalizer. The filter tap coefficients are analyzed and the sampling phase Is adjusted such 
that a cost function associated with the performance of the equalizer structure Is substantially optimized. Unfortunately, 
for very distorted channels, feed forward equalizers and decision feedback equalizers fall due to their high BER. Rather 
Vitorbl equalizers should be used. Generally, the methods based on symbol rate sampling take a long lime for syn- 
chronization. 

[0015] Oversampling methods perform sampling of the received signal at a sampling frequency higher than symbol 
frequency. 

[0016] US 5,291,523 titled * Viterbi Receiver with Improved Timing Means* teaches to use two Viterbi detectors and 
sample the received signal twice in each symbol period. 

[0017] In accordance with the disclosure of EP 1 009 125 A2 titled "Viterbi Receiver With Improved Timing Means", 
synchronization with the received high-speed serial digital signal can be establish ed at high-speed, and synchronization 
can be adjusted towards the zero transition point of the received signal. Fig. 3 of EP 1 009 125 A2 Illustrates four-fold 
over-sampling. 

[0018] A received symbol stream Is M-fold and 12-fold oversampled and correlated with a known synchronization 
word in US 5,533,066 ("Apparatus and Method for Estimating Maximum Likelihood Sequence Using Optimum Sampling 
Phase") and WO 02/076010 A2 ("Symbol Recovery From an Oversampled Hard-Decision Binary Stream"), respec- 
tively. However, usage of a synchronization word Is undesirable in many applications. Moreover, the main weakness 
of techniques using massive oversampling Is the required high processing speed. N-fold oversampling necessitates 
processing at a frequency N times higher than symbol frequency. This Is not feasible in high-speed systems such as 
10 GHz or 40 GHz optical systems. 

[0019] In WO 02/30035 A1 the Inventors describe symbol timing in a system which does not provide a carrier cor- 
responding to a symbol frequency. By coilecting a histogram of samples for a predetermined number of symbol times 
symbol edges and a maximum eye opening are determined. Specifically an average, weighted average, or other method 
Is applied to determine an average timing for maximum eye opening for each symbol time, Eight-fold over-sampling 
is employed. 

[0020] It is the object of this invention to provide an improved method and an Improved sampling phase adjustment 
circuit. 

[0021] This object Is achieved by the subject matter of the Independent claims. 

[0022] Preferred embodiments of the invention are the subject matter of the dependent claims. 

[0023] The advantage of choosing a sampling phase so that a population difference parameter Is maximized is that 

this leads to a, at least nearly, optimum, i.e. to a minimized bit error rate. Following the de-facto conventions In the 

literature, the terms "optimal" or "optimized" are used In a somewhat loose sense. What is meant is that a solution of 

minimized BER is sought within some practical framework or solution space, not excluding the case that in a slightly 

modified framework even lower BER might be achieved. 

[0024] Provided that the first and second digital values may assume N different digital values, a number of N counters 
In each of the first and second plurality of counters provides the most accurate approximation of the population differ- 
ence parameter defined In equation (1). 
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[0025} Asimple and robust algorithm for maximizing the population difference parameter and thereby optimizing the 
sampling phase in order to arrive at a minimum BER Is provided. To this end the sampling phase is increased, if the 
population difference parameter increases with the sampling phase and decreased otherwise. More specifically, the 
sampling phase Is increased if asecond population difference parameter following afirst population difference param- 
eter is bigger than the first population difference parameter and the sampling phase has been Increased between the 
first and second population difference parameter. The sampling phase Is also Increased if the second population dif- 
ference parameter Is smallerthanthefirstpopulation difference parameter and the sampling phase has been decreased 
between the first and second population difference parameter. 

[0026J Using dedicated circuits like counters, absolute difference circuits, de-multiplexer and first and second distri- 
bution circuits provide for a high operation speed. 

[0027] Since the counters count M quantized samples the necessary operation speed after the counters drops ap- 
proximately by a factor M wherein M may be chosen to be two thousand. As a consequence general purpose logic 
circuits like a microprocessor or CPU may be used rather than dedicated logic In order to save development and 
production costs for low volume batches. 

[0028] In order to avoid a delay when a microprocessor or CPU sequentially reads out the counter values holding 
circuits may be provided between the counters and the microprocessors for temporarily storing the counter values after 
M samples have been processed, So the holding circuits constitute an Interface which Insures that ail quantizedsamples 
are taken into consideration for the calculation of the population difference parameter i.e. no quantized samples are 
Ignored during an idle period after M quantized samples have been processed. 

[0029] This invention can be more fully understood from the following detailed description when taken In conjunction 
with the accompanying drawings, in which: 

Fig. 1a shows one embodiment of an Inventive sampling phase adjustment circuit; 
Fig.1 b shows a second embodiment of an inventive sampling phase adjustment circuit; 

Fig. 2 shows some relevant units of a digital receiver wh Ich may bo used for high speed optical data transmission; 
and 

Fig, 3 shows the dependence of the bit error rate and a population difference parameter on the sampling phase. 

Abbreviations 

[0030] 

ADC; Analog-to-dlgital converter 
AGC: Automatic gain control 
BER; Bit error rate 
CPU; Central processing unit 
CR: clock recovery 

ECC: Error correcting code, error correction circuit 

EQ: digital equalizer 

FSE: Fractionally spaced equalizer 

ISI: intersymbol Interference 

PS: phase shift 

SPA: sampling phase adjustment 

XOR: exclusive-OR 

Mathematical Symbols 

[0031] 

b(t): timing signal 

cip k (l): counter state of counter i for phase <p k after M symbols 

c(t): control signal 

c!k: ciock 

n: resolution in bits of the quantizer 

N: number of counters In each group, number of quantization levels 
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M: 



number of symbols counted 



pdf: 


probability density function 


pd: 


population difference parameter 


pdm: 


modified population difference parameter 


Pd q r 


quantized population difference parameter 


pdm q : 


quantized modified population difference parameter 


T: 


symbol period 


A: 


sampling phase adjustment step 


<p: 


phase 




sampling phases 



[0032] While the present invention Is described with reference to the embodiments as illustrated In the following 
detailed description as well as In the drawings, it should be understood that the following detailed description as well 
as the drawings are not intended to limit the present Invention to the particular illustrative embodiments disclosed, but 
ratherthe described illustrative embodiments merely exemplify the various aspects of the present Invention, the scope 
of which Is defined by the appended claims. 

Sampling Phase Adjustment Method 

[0033] The sampling phase adjustment (SPA) method exploits sample statistics in order to adjust the sampling In- 
stants to their correct value. The sample statistics are obtained from the waveform of a plurality of symbols or pulses 
at different sampling instants. The SPA method is performed within a receiver of digitally transmitted signals. The 
relevant parts of the receiver are presented in Fig. 2. The received electrical signal is designated s(t). It may be the 
output signal of a preamplifled antenna signal In the case of mobile communications or the output of a optical-Jo- 
electrlcal Interface of an optical receiver. Due to varying power level of the received signal s{t), an Automatic Gain 
Control (AGC) circuit 1 is included before the Analog-to-Digital Converter (ADC) 4 in order to generate a gain controlled 
signal gs(t). The AGC circuit 1 ensures a proper analog-to-dlgital conversion of the analog signals. ADC 4 oulpuls 
quantized samples qs. A digital equalizer (FQ) 5 and an error correcting circuit (ECC) 6 usually follow the ADC. 
[0034] A clock recovery (CR) circuit 2 extracts frequency and phase and generates a local sampling clock at equi- 
distant points in tlmo, with some fixed phase relation to the transmitted symbol stream. As explained above (cf. 
Buchwald95), Cfl circuit 2 may be connected to the output of the AGC 1 In order to recover the clock on the bases of 
the gain controlled signal gs{l). Alternatively or In addition the CR circuit 2 may be connected via CR logic circuit 8 to 
the output of ADC 4 In order to recover the clock on the basis of tho quantized samples (cf. US 5,048,060). 
[0035] A Sampling Phase Adjustment (SPA) circuit 3 adds a phase delay to the recovered sampling clock, In order 
to optimize receiverperformance. The phase delay Is determined by SPA logic circuit 7 based on the quantized samples 
qs. Two embodiments of SPA logic circuit 7 are shown in Fig, 1a and 1b. 

[0036] In order to reduce the bit error rate (BER) at tho output of equalizers, which is implemented as Viterbi equalizer 
in one embodiment, the signal is oversampled at twice the symbol rate. 

[0037] Samples representing Ihe transmitted symbols are processed to get the correct sampling phase. The two 
samples representing a transmitted symbol are sampled at phase tp k , k ={1 ,2}, 0 < tp„ S 2n. To find Ihe best sampling 
phase a population difference parameter pd defined by equation (1) may be utilized: 



pdf denotes the amplitude probability density function at either oversampllng phase, pd Is called a population difference 
parameter, because it describes the overall difference in population (density) per possible amplitude level. 
[0038] There is no practical way for calculating the values of these probability density functions in real systems. 
Consequently, these distributions are approximated by their discrete counterparts, i.e. by amplitude histograms ob- 
tained from processing quantized samples. 

[0039] The ADC 4 has N quantization levels. In an embodiment N Is a power of 2: N=2 n . in this embodiment n bits 
are used to represent a sample value. The N different amplitudes are mapped into digital words denoted by a,, i = 
{1,2„..,N} , and a f =i-1. For a predetermined duration the timing recovery circuit counts the occurrences of each possible 
word at the two different sampling instants <p k , k={1 ,2} by using 2N counters. By applying this counting, two histograms 
approximating the probability density functions are obtained. The counter states after M symbols are denoted ctp k (i). 




(1) 
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A quantized population difference parameter pdq is defined in equation (2): 




(2) 



£ % (/)=Mfork={1,2} 



(3) 



[0040] The sum of ail counter states Ccp^l ) after M symbols Is M as indicated by equation (3). For reasons of practical 
Implementation, we Introduce a quantized modified population difference parameter pdmq: 



Implementation 

[0041] An embodiment of the SPA logic circuit 7 for calculating pdmq and timing adjustment is shown in Fig. 1a. It 
comprises a de-multiplexer 11, distribution circuits 12 and 13, counters 14, 15, 16, 17, 18, 19, 20 and 21, absolute 
difference circuits 22, 23, 24 and 25, adder 26, a holding circuit 27, a comparator 28, an excluslve-OR (XOR) circuit 
29 and a phase calculation circuit 31 , The quantized samples qs are first demultiplexed by demultiplexer 11, The 
quantized samples qs, which have been sampled at phase (f>, are provided to distribution circuit 12 by demultiplexer 
1 1 . On the other hand quantized samples qs 2 sampled at phase tp 2 areprovlded to distribution circuit 1 3 by demultiplexer 
1 1 . The distribution circuits compare each quantized sample to a plurality of digital values the quantized samples may 
assume, To each of the digital values a counter is associated, The distribution circuit generates a clock pulse for every 
quantized sample it receives In order to Increment the counter to which Ihe digital value of the quantized sample is 
associated with, 'thereafter the 2N counters 1 4 to 21 count all events during a predetermined time period MT, where 
T is the symbol period. The absolute difference circuits 22 to 25 and the adder 26 Implement equation (4) i. e. calculate 
the pdm q (t) parameter value at time t. 

[0042] In another embodiment the ADC 4 may already have two outputs, a first output for the quantized samples 
qs, which have been sampled at phases, and a second output fortho quantizedsamplesqsj which have been sampled 
at phase ife, This embodiment does not comprise demultiplexer 11 , Rather the first and second outputs are directly 
connected to the Inputs of distribution circuits 12 and 13, respectively. 

[0043] Afurtherembodimentcomprisesa1:m-typedemuitiplexer 11, m distribution circuits, m counter groups, each 
comprising N counters and N/2 absolute difference circuits. Consequently the circuitry between demultiplexer 11 and 
adder 26 is m-fold parallel. Provided that 2-fold oversampling is employed In this embodiment, the m-fold parallelism 
decreases the necessary processing speed of the distribution circuits, counters and absolute value circuits by a faclor 
of nV2 compared to the embodiment shown In Fig. 1 a. 

[0044] Holding circuit 27 stores the previously calculated pdm q {t-MT) parameter value which has been calculated at 
time t-MT. Comparator 28 compares the newly calculated pdm q {,t) parameter value with the previously calculated pdm q 
(t-MT) parameter value in order to output control signal c(t) e {-1,1}. Control signal c(t) Is input Into XOR circuit 29. 
XOR circuit 29 generates a timing signal b(t)e {-1 ,+1} to be used for sampling phase adjustment. Timing signal b(t) Is 
feed back to the other input of XOR circuit 29 via holding circuit 30 in order to delay timing signal b(t) by MT. Conse- 
quently the XOR circuit 29 XORs c(t) and b(t-MT). Tabie 1 depicts its functionality. 




(4) 



Table 1: 



functionality of XOR circuit 29 



c{t) b(t-MT) 



b(t) 
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Table 1: (continued) 



functionality of XOR circuit 29 


c(t) 


b(t-MT) 


b(t) 


- 1 


1 


■1 


1 


-1 


-1 


1 


1 


1 



[0045] Practically, the sampling phase adjustment circuit 3 works in a dithering mode and during an adjustment time 
this circuit will tend to find the desired sampling phase, that is, ifiesampllng phase maximizing the population difference 
parameter pdmq. 

[0046] The new value of a sampling phase <p(t) is obtained by using equation (5): 



cp(f)^<r-/wr) + b(f)A (5) 



[Q047J in equation (5) (p(t-MT) represents the previous or old sampling phase and A represents a sampling phase 
adjuslment step. Equation (5) is implemented In phase calculation circuit 31 . Sampling phase p(t) Is provided to SPA 
circuit 3 which updates the phase shift <p(t) In compliance with equation (5) after a time of MT and shifts the clock elk 
(t) provided by clock recovery circuit 2 by (<p(t)/2jt)T in order to generate a shifted clock clk(t+(<p(t)/2n)T). After a few 
phase adjustments, the sampling phase providing the maximum value of the pdm q parameter will be reached, and 
subsequently be tracked. 

[0048] In an alternate embodiment XOR circuit 29 may be a conventional logic circuit inputting either logic 1s or 0s 
and outputting Os or 1s. In this embodiment equation 5 must be modified in that previous sampling phase (ji(t-MT) is 
increased by A if b(t) Is one and old sampling phase is decreased by A if b(t) Is equivalent to 0 In phase calculation 
circuit 31. 

[0049] Fig. 1 b shows a second embodiment of SPA logic circuit 7. in contrast to the first embodiment shown !n Fig, 
1 a the elements 22 to 31 are replaced by holding circuits 41 , 42, 43, 44, 45, 46, 47 and 48, and microprocessor 49. 
Upon processing M symbols by de-multiplexer 11 , distribution circuits 12 and 13 and counters 14 to 21 , the counter 
values are stored In holding circuits 41 to 48. Then the counters are reset and are ready for processing the next 2M 
samples of the next M symbols. Microprocessor 49 then reads out the values of the holding circuits 41 to 48 and 
calculates a population difference parameter pd q or pdm q in compliance with equations 2 or 4, respectively. Then the 
microprocessor performs the functionality of holding circuit 27, comparator 28, XOR circuit 29, holding circuit 30 and 
phase calculation circuit 31 as described in connection with Fig. la. 

[0050] As wilt be explained later In connection with the simulation result a typical value for M Is 2000. As a conse- 
quence the required operating frequency of the multiplexer 11, distribution circuits 12 and 13 and counlers 14 to 21 
must have a substantially higher operation speed than microprocessor 49. Although microprocessor 49 has to perform 
more than one operation In order to calculate equation (2) or (4) and perform the maximum search a substantial dif- 
ference In operation speed of roughly a factor 20 remains for M equivalent to 2000 provided that the processing of the 
counter values takes 1 00 microprocessor clock cycles, So even for high speed applications a general-purpose micro- 
processor 49 having a moderate clock frequency can be used, The application of standard microprocessors reduces 
development costs and receiver costs especially for small batch sizes, 

[0051] The holding circuits 41 to 48 may be omitted If microprocessor 49 is quick enough to enter all counter values 
within a short time or a certain number of samples are discarded in an idle period between the blocks of M samples 
which are counted. In the Idle period the microprocessor 49 may read out the counter values. 
[0052] In another embodiment the counter values are read out by CPU one after another in a distance smaller than 
MT/(2N). Then each counter value is reset. In this embodiment the periods of time each counter evaluates are shifted 
with respect to each other. Since it is assumed that the symbols do not appear In a certain pattern but rather are 
statistically Independent this shift may be acceptable. 

[0053] In a third embodiment which is not shown In a figure the quantized samples may be directly input into a 
microprocessor which performs the function of elements 11 to 31 shown In Fig. 1 a. Such a circuit structure can be 
applied only to moderate transfer rates due to operation speed limitation of the microprocessor 49. In this embodiment 
counters 14 to 21 may be Implemented slmpty by registers which are increased by an adder logic. 
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Simulation 

[0054] To show the power of the Inventive method art optical channel with a high amount of chromatic dispersion, 
polarization mode dispersion, and noise was simulated. The sampling resolution was 1/32 equivalent to 32-fold over- 

5 sampling for simulation purposes In orderto arrive at a quasi-continuous signal. After ADC 4, a fractionally spaced ML 
equalizer using two samples in a bit slot for a branch metric calculation was applied. The BER obtained at the output 
of the fractionally spaced ML equalizer is shown In Fig. 3, where the bit error rate shown at time instant t is obtained 
by using samples at point t and t+T/2. In Fig. 3 the pdm q parameter Is represented, normalized by Its maximum value. 
It can be seen that sampling at time Instants (t f , t 2 ) produces a nearly optimum BER performance of 4.10" 4 . In fact, 

10 the true minimum Is 4/32T earlier; however, for practical purposes, this deviation from true minimum is not significant. 
In this connection reference Is made to Fig. 2, specifically ECC 6 which can operate at BEfls of 1 0' 3 and reduces these 
BERs to 10" 9 to 1 0- 15 at its output. Nearly optimum BER is also obtained for the case of sampling at time Instants (tg, 
t 3 ). Sampling at time instants ft,', t 2 ") leads to satisfactory BER, which, however, is still significantly worse than the 
BER produced when sampling at time instants (t 1t t 2 ). 

is [0055] in one bit slot there are two local minima and two tocal maxima of the pdmq parameter. By selecting either 
of these maxima the equalizer 6 and ECC 6 will have almost minimum BER at its respective output. 
[0056] Auxiliary simulation results indicate that aboul 2000 symbols are enough to get a sufficiently accurate value 
of the pdmq parameter for use in sampling phase adjustment. 

[0057] Further modifications and variations of the present Invention will be apparent to those skilled In the art in view 
so of this description. Accordingly, this description Is to be construed as Illustrative only and Is for the purpose of teaching 
those skilled in the art the general manner of carrying out the present invention. It is to be understood that the forms 
of the invention shown and described herein are to be taken as the presently preferred embodiments. 



SB Claims 

1. A self-timing method for phase adjustment of a sampling phase (CLK) to an analog signal (s(t)) comprising a 
stream of symbols, the method comprising: 

30 digitizing (4) Ihe analog signal at a first phase with respect to said symbols In order to obtain a first quantized 

sample (qs,) having a first digital value out o! a first plurality of digital values; 

digitizing (4) the analog signal at a second phase with respect to said symbols in order to obtain a second 
quantized sample (qs 2 ) having a second digital value out of a second plurality of digital values; and 

35 

repeating said digitizing at said first and second phases for different symbols comprised in said analog signal; 
characterized by: 

40 Increasing a first counter out of a first plurality of counters (14, 16, 1 6, 17); said first counter being associated 

with said first digital value; 

Increasing a second counter out of a second plurality of counters (18, 19, 20, 21); said second counter being 
associated with said second digital value; 

45 

said repeating further comprises repeating said Increasing of said first and second counters; 

adjusting (31) said sampling phase based on the values of said counters of said first and second plurality of 
counters. 

50 

2. The method of claim 1 , characterized in that a population difference parameter (pdm q ) is calculated based on 
the values of the counters of said first and second plurality of counters and said sampling phase being adjusted 
that s aid population difference param eter is maximized. 

55 3. The method of one of the claims above, characterized in that said first and second pluralities of counters comprise 
N counters each. 



4. The method of claim 3, characterized In that each counter in said first plurality of counters Is associated with a 
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counter In said second plurality of counters; a population difference parameter is the sum (26) of the absolute 
values of the differences (22, 23, 24, 25) between the values of said counters of said first plurality of counters and 
the values of the corresponding counters in said second plurality of counters. 

The method of one of the claims above, characterized In that said first and second plurality of digital values 
comprise N different values and said first and second pluralities of counters comprise N counters each. 

The method of claim 2 or claims 3 to 5, as far as they refer to c!alm 2, characterized further by: 

repeatedly calculating said population difference parameter, 

comparing (28) a first population difference parameter (pdm q (t-MT)) to asecond population difference param- 
eter (pdm q (t)) following said first population difference parameter; 

increasing (31) said sampling phase ((p(t)), If said second population difference parameter (pdm q (t)) is bigger 
(28) than said first population difference parameter (pdm q (t-MT)) and said sampling phase has been increased 
(b(«-MT)) after the calculation of said first population difference parameter or if said second population differ- 
ence parameter (pdm q (t)) Is smaller than said first population difference parameter (pdm q (t-MT) and said sam- 
pling phase has been decreased (b(t-MT)) after the calculation of said first population difference parameter 
(pdm q (t-MT)); and 

decreasing said sampling phase (<p(t)) if said second population difference parameter (pdm q (t)) is smaller (28) 
than said first population difference parameter (pdm q (t-MT)) and said sampling phase has been Increased 
after the calculation of said first population difference parameter (b(t-MT)) or If said second population differ- 
ence parameter Is bigger (28) than said first population difference parameter (pdm q (l-MT)) and said sampling 
phase has been decreased after the calculation of said first population difference parameter (pdm q (t-MT)). 

A digitizing self-timing circuit comprising: 

an analog- to-digital converter (4) receiving an analogue signal (gs(t)) comprising a stream of symbols; said 
analog-lo-dlgilal counter (4) performing two analog-to-dlgltal conversions per symbol period (T), thereby pro- 
ducing first quantized samples (qsi) at a first phase and second quantized samples (qs ? ) at a second phase; 
said first quantized samples have digital values out of a first plurality of digital values and said second quantized 
samples having digital values out of a second plurality of digital values; 

characterized by: 

a first plurality of counters (14, 15, 16, 17) each counter of said first plurality being associated with a first digital 
value and being Incremented if a first quantized sample (qS,) Is equivalent to said first digital value; 

a second plurality of counters (18, 19, 20, 21) each counter of said second plurality being associated with a 
second digital value and being incremented if a second quantized sample (qs 2 ) Is equivalent to said second 
digital value; and 

a sampling phase adjustment circuit (3) for controlling the sampling times of said analog-to-digltal counter (4) 
based on the values of the counters of said first and second plurality of counters. 

The circuit of claim 7, further comprising: 

circuitry (22, 23, 24, 25, 26) connected to said counters of s aid first and second plurality of counters calculating 
a population difference parameter (pdm q ) from the values of the counters of said first and second plurality of 
counters and a maximizing circuit (27, 2B, 29, 30, 31) being connected to said circuitry (22, 23, 24, 25, 26) 
and said sampling phase adjustment circuit (3) for controlling said sampling phase adjustment circuit in order 
to maximize said population dilference parameter (pdm,,). 

The circuit of claim 8, each counter of said first plurality of counters being associated to a different one of said 
second plurality of counters in a one to one relationship, said circuitry comprising: 
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a plurality of absolute difference circuits (22, 23,24, 25) each absolute difference circuit being connected to a 
pair of associated counters and calculating the absolute difference between the counter values of said pair of 
associated counters; and 

an adder (26) connected to the outputs of said plurality of absolute difference circuits and said maximizing 
circuit (27, 28, 29, 30, 31) for summing the absolute differences and providing the sum to said maximizing 
circuit (27, 28, 29, 30, 31 ). 

10. The circuit claim 7 or 8, said maximizing circuit further comprising: 

a comparator circuit (27, 28) for receiving the output of said circuitry (22, 23, 24, 25, 26) and comparing a first 
population difference parameter (pdm q (1-MT)) to a second population difference parameter (pdm q (t)) following 
said first population difference parameter (pdm q (t-MT); 

a phase calculation circuit (29, 30, 31) connected to the output of said comparator circuit (27, 28) and to said 
sampling phase adjustment circuit (3) for outputtlng a phase shift to said sampling phase adjustment circuit 
(3); said phase calculation circuit (29, 30, 31) 

Increasing said sampling phase ((p(t)) if said second population difference parameter (pdm q (t)) Is bigger 
(28) than said first population difference parameter (pdm q {t-MT)) and said sampling phase has been increased 
(b(t-MT)) after the calculation of said first population difference parameter or if said second population differ- 
ence parameter (pdm q (t)) is smaller than said first population difference parameter (pdm q (t-MT)) and said 
sampling phase has been decreased (b(t-MT)) after the calculation of said first population difference parameter 
(pdm q (t-MT)}; and 

decreaslngsaldphasecalculatlon circuit (29, 30, 31) decreasing saldsBmpllngphase(ij>(t))tfsaldsecond 
population difference parameter (pdm q (t)) is smaller (28) than said first population difference parameter (pdm q 
(t-MT)) and said sampling phase has been increased (b(t-MT)) after the calculation of said first population 
difference parameter or if said second population difference parameter (pdm q (t)) Is bigger than said first pop- 
ulation difference parameter (pdm q (t-MT)) and said sampling phase has been decreased (b(t-MT)) after the 
calculation of said first population difference parameter (pdm q (t-MT)). 

11. Theclrcuitofoneof claims 7 to 9, f u rther comprising : 

a dc-mulliplcxor (11 ) connected to said analcg-to-rilgltal converter (4) for outputtlng said first quantized sam- 
ples (qS() at a first output and said second quantized samples (qs 2 ) at a second output; 

a first distribution circuit (12) connected to said first output of said de-multiplexer (11) and each of said first 
plurality of counters; said first distribution circuit comparing a first quantized sample to a plurality of digital 
values to each of which a counter of said first plurality of counters (1 4, 1 5, 1 6, 1 7) Is associated and generating 
a clock pulse to Increment the counter to which the digital value is associated which is equivalent to said first 
quantized sample; and 

a second distribution circuit (13) connected to said second output of said de-multiplexer(11) and each of said 
first plurality of counters; said second distribution circuit comparing a second quantized sample to a plurality 
of digital values to each of which a counter of said second plurality of counters (1 8, 1 9, 20, 21 ) is associated 
and generating a clock pulse to Increment the counter to which the digital value is associated which is equiv- 
alent to said second quantized sample. 
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(54) Method and circuit for controlling amplification 



(57) This invention relates to a method for control- 
ling an amplification within a receiver. The gain of a var- 
iable gain amplifier which amplifies an Input slgnai to 
obtain a gain-controlled signal is set. Said gain control- 
led signal Is anaiog-to-dlgltal converted. Thereby digital 
words are generated. Each digital word has a value out 
of a plurality of possible digital values. Digital words hav- 
ing a value within a first and a second subset of the pos- 



sible digital values are counted In order to generate a 
first and a second counter value. The gain is set in ac- 
cordance with the first and second counter values In a 
fash Ion that all counter values are as equal as possible. 
Furthermore, this invention is related to a circuit to be 
used within a recoivor for adjusting tho gain of a variable 
gain amplifier in order to ensure a proper analog-lo-dlg- 
ital conversion. 
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Description 

[0001] The present invention relates to a method and a circuit according to the preamble parts of claims 1 and 6, 
respectively, in particular the invention relates to the automatic gain control In digital receivers e.g. for optical data 
5 transmission. 

[0002] Such a method and circuit are described In US 3,931 ,584 titled "AUTOMATIC GAIN CONTROL". 
[0003] Communication systems are designed to provide reliable transmission of information, whereby they should 
be able to combat various channel imperfections, in digital communication systems, the most Important parameterfor 
the signal quality measurement surely is BER (bit error rate). The signal at the receiver end contains the data and dock 

10 Information to be extracted. The amplitude of the received signal is attenuated, whereby the attenuation depends on 
length and type of the transmission channel. Due to varying power level of both wired and wireless digital transmitted 
signals, an automatic gain control circuit (AGC) Is customary placed before a CDR (clock and data recovery) circuit, 
which comprises a clock recovery circuit (CR), an anaiog-to-dlgital converter (ADC) and a sampling circuit, which 
receives the recovered clock from the ciock recovery circuit, 

's [0004] Automatic gain control is widely employed to adjust the amplification of an Input signal In such away that an 
output signal has a substantially constant signal level In terms of root mean square, mean or peak values, Most digital 
communication systems must accommodate received signal levels varying within a range of some tens of dB. However, 
the power level variation of the AGC output signal should be as low as possible and at (east within a range of less than 
one dB. 

so [0005] in US 3,931 ,584 Inventors present an AGC separated In two circuits one for coarse and the other for fine gain 
control. In order to adjust the amplification of the tine AGC, the sign and most significant bit in a digital word that an 
ADC circuit produces are monitored. An exclusfve-OR (XOR) operation is performed on the sign and most significant 
bit In order to detect zero-one and one-zero combinations which Indicate high positive and low negative values In a 
two's compliment system, Then the zero-one and one-zero combinations are counted within a period of time. According 

ss to this document It will be apparent to those skilled In the art that any combination of the digits can be counted to control 
the amplification of the AGC. A comparator compares the count with a threshold value. If the count Is less than the 
threshold value the fino AGC increases the power level and vice versa. Thereby amplification of the input signal is set 
to a level that the highest and lowest digital words approach the maximum and minimum, respectively, level expressible 
by the digital words. So the quantization noise of the converter is reduced by maximizing the number of digital words 

so that are available to express the signal amplitude, 

[0006] This method has the shortcomings that the th rcshold value depends on the communication channel conditions 
such as noise power, noise distribution, dispersion etc. In case of very distorted channels like wavelength division 
multiplexing (WDM) or wireless channels the conventional methods might provide sub-optimal or even very bad per- 
formance. In time-variant channels, the measurement of channel performance and selection of new thresholds for 

35 obtaining optimal decoding performance increase the AGCcircult complexity. These measurements are hardly feasible 
for very distorted channels. The receiver developed for such channels have to utilize either a channel equalizer or 
strong ECC (error correcting code) or a combination thereof. 

[0007] It is the object of this Invention to present a method and circuit for controlling amplification within a receiver 

In order to arrive at an optimized BER. 
40 [0008] This object Is achieved by the subject matter of the Independent claims. 

[0009] Preferred embodiments of the invention are the subject matters of the dependent claims. 

[0010] Following the de-facto conventions In the literature the terms "optimal" or "optimized" are used in a somewhat 

loose sense. What is meant is that a solution of minimized BER Is sought within some practical framework or solution 

space, not excluding the case that In a slightly modified framework even lower BER might be achieved. 
4S [001 1] Grouping more than one possible digital value Into subsets advantageously reduces the number of counters 

required and consequently reduces circuit complexity. 

[0012] Minimizing a uniformity parameter provides an efficient method for adjusting the amplification to achieve a 
close to minimum BER even for highly distorted, time-variant signals. This technique can be applied to non-Gaussian 
time-variant channels in the presence of linear and nonlinear distortion. 
so [001 3] Defining a first, inner subset and a second, outer subset and comparing the counter values of the Inner and 
outer subsets avoid a minimum search. The amplification Is increased If the Inner subset Is greater than the counter 
value of the outer subset, 

[001 4] The difference of the counter value of the inner subset minus the counter value of the outer subset constitutes 
a control parameter, which can be advantageously used together with popular control algorithms like a proportional- 
55 plus-lntegral-ptus-derivative controller (PID controller) or fuzzy controllers. 

[001 5] In the following preferred embodiments of this invention are described ref erringto the accompanying drawings. 
[0016] Figure 1 shows a block diagram of the relevant units of a receiver according to this invention, 
[0017] Figure 2a shows a detailed block diagram of a first embodiment of the fine logic unit. 
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[0018] Figure 2b shows a detailed block diagram of a second embodiment of the fine logic u nit. 

[0019] Figure 3 shows a statistic of digitized sample values. 

[0020] Figure 4 Illustrates counter values for the statistics of figure 3 for a first subset. 

[0021] Figure 5 Illustrates counter values for the statistics of figure 3 for a second subset. 

s [0022] Figure 6 shows simulation results. 

Abbreviations 



[0023] 
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ADC: 


analog-to-digital converter 


AGC: 


automatic gain control 


BER; 


Bit error rate 


CDR: 


clock and data recovery 


ts CR: 


clock recovery 


ISI: 


Intersymbol Interference 


WDM: 


wavelength division 




multiplexing 


XOR: 


excluslve-OR 


so ECC: 


Error correcting code, error correction circuit 


EQ: 


digital equalizer 


MLSE 


Maximum-likelihood sequence estimator 


VGA: 


variable gain amplifier 


25 Mathematical Symbols 


[0024] 




a,: 


digital word 




quantized sample at time j 




count of subset S q 


symbol rate 


fg: 


fine gain control signal 



I: digital word Index, 1 e i <; 2° 
ss j : time Index, 1 s M 
M: number of samples 
k q : expected number of samples In subset S q 
I: number of subsets 
m e : coarse gain parameter 
40 m u : uniformity parameter 
n: ADC resolution in bit 
p: index for uniformity parameters m u 
q: subset Index, 1 < q < I 
S: complete set 
45 s 0 : coarse gain subset 
S q : subset 
t: time 
T: period sample 
T t : Idle period 

so 

[0025] While the present invention is described with reference to the embodiments as Illustrated in the following 
detailed description as well as in the drawings, it should be understood that the following detailed description as well 
as the drawings are not intended to limit the present invention to the particular Illustrative embodiments disclosed, but 
rather the described illustrative embodiments merely exemplify the various aspects of the present Invention, the scope 
ss of which is defined by the appended claims. 

[0026] Several counters are introduced, some used for the coarse AGC unit and some others used for fine AGC 
adjustment, as described in the following. 

[0027] Since the gains of the coarse and fine AGC units are controlled by other elements It may be more appropriate 



3 



EP 1 473 831 A1 



to designate these units as variable gain amplifiers (VGA). AGC and VGA will be used synonymously in the following. 
Implementation 

5 [0028J An excellent AGC capability can be realized by using several AGC units. Typically a coarse AGC unit has a 
wider dynamic rage than a fine AGC unit, which provides for a more accurate gain control, in one embodiment, the 
coarse AGC unit Is controlled In accordance with conventional practice. A peak detector detects either only positive 
peaks or positive and negative peaks In the output signal of the coarse AGC unit and adjusts its amplification accord- 
ingly. In particular it is detected when that the output signal of the coarse AGC unit exceeds a positive threshold or 

to exceeds the positive threshold and falls below a negative threshold for a predetermined portion of time. 

[0029] In another embodiment, which is described In connection with Fig. 1 , the coarse AGC unit can also be con- 
trolled by signal generated in a logic unit or a DSP (digital signal processing) circuit. At start-up, the coarse AGC unit 
can quickly bring a total gain into working regime and after that the fine AGC would adjust finally the gain providing 
the minimum BER at the output of the digital equalizer or other decoding circuit. 

15 [0030] More specifically, after power-up the gain of the fine AGC is set to some relaxed initial value, which Is an 
Intermediate gain value that Is neither ciose to the maximum nor to the minimum gain of the fine AGC unit in order to 
allow positive and negative gain adjustments. Then, In an coarse adjustment phase only the gain of the coarse AGC 
unit Is adjusted to bring the Input signal of the following ADC into She working regime i. e. a coarse gain parameter m 0 
remains above some lower threshold. When the Input signal has reached the working regime, the gain of the coarse 

so AGC unit remains fixed whereas the gain of the fine AGC unit optimizes the level of the Input signal to arrive at a 
minimum BER during normal operation, The receiver may reenter Into the coarse adjustment phase in order to re- 
adjust the coarse AGC unit, when the gain of the fine AGC unit reaches either its maximum or minimum. During the 
coarse adjustment phase the gain of the coarse AGC unit Is adjusted In such a way that a coarse gain parameter m 0 
always remains above some lower threshold. The most difficult situation expected In an optically amplified DWDM 

25 system is when many channels are switched-away or are added at the same time In an optical add/drop multiplexer 
(ADM) which causes large power transients in optical amplifiers. When re-entering Into the coarse adjustment phase, 
the gain of the fine AGC is again set to its relaxed Initial value and the gain of the coarse AGC unit is adjusted to bring 
the total gain into the working regime as explained above, 

[0031] Figure 1 shows a block diagram of the relevant units of a receiver according Jo this Invention. It comprises a 
so coarse VGA unit 2, a fine VGA unit 3, a clock recovery (CR) 4, a sampler 5, an ADC 6, an digital equalizer 7 and a fine 
logic unit 8 and a coarse logic unit 9, The Input signal 10 that Is input into the coarse VGA unit 2 represents the 
transmitted signal as distorted by the channel at the receiver end. It may be the output signal of an antenna In the case 
of radio communication or the output signal from an opt!cal-to-electrical interface in the case of optical data transmis- 
sion. 

35 [0032] The signal at the input of the clock recovery circuit 4 must be strong enough. However, it Is even more critical 
to provide sampler 5 and ADC 6 with a signal having optimal strength In order to arrive at a minimum BER, To this end 
coarse and fine VGA 2 and 3 amplify the Input signal and generate a gain controlled signal 12 which is Input into CR 
circuit 4 and sampler 5, CR circuit 4 recovers the symbol clock and provides control signals on lines 13 and 11 to 
sampler 5 and ADC 6. These control signals already have an appropriate phase relation to the gain-controlled signal 

40 12. Such ciock recovery circuits are described In "Integrated Fiber-Optlc Receivers" by A. Buchwald, K. W. Martin, 
Kluwer Academic Publishers, 1995. It is noted that ADC 6 and sampler 5 could be reversed in sequence. 
[0033] The sampler 5 essentially comprises a samplc-and-hold circuit, which samples the gain-controlled signal 12 
at appropriate points In time and outputs a constant signal on lino 1 4 to the ADC until the next sampling is performed. 
The ADC provides Its output In form of digital words on line 15 to digital equalizer 7 and fine and coarse logic unit 8 

& and 9. The digital equalizer 7 preferably combats ISI and consequently provides equalizer functionality. The fine ioglc 
unit 8 essentially calculates statistics of the digital words provided by the ADC 6 and generates a fine gain control 
signal on tine 1 8 to control the gain of fine VGA 3. The coarse logic unit 9 operates in a similar fashion and generates 
a coarse gain control signal on line 19 to conlrol the gain of coarse VGA 2. The operation of the coarse and fine logic 
unit 8 and 9 will be explained in more detail below. Two embodiments of fine ioglc unit 8 are shown In figures 2a and 2b. 

so [0034] In another embodiment sampler 5 and ADC 6 may perform over-sampling. As a consequence, more than one 
digital words are provided on line 15 for each symbol comprised In the input signal 10, In this case, digital equalizer 7 
may reduce the over-sampled digital words to one digital word per symbol and provide this data on line 20 to logic 
units 8 and 9 or alternatively the over-sampled digital words.may be processed directly by the logic units, 
[0035] ADC 6 has 2" quantization levels, where each sample Is represented by n bit. The 2 n different quantization 

65 levels are mapped Into a set S of digital words, denoted by S={a|}, 1=1 , 2, ... , 2". Let the subscript i of the digital word 
a, represent the amplitude level of the quantized sample In value order, I.e. a 1 <a 2 <...<a 2 n . 

[0036] In a particular interval of time MTthe DSP counts the occurrences of each possible word, by using 2 n counters. 
Here M Is the number of all words In this Interval and T Is the symbol period being the reciprocal of the symbol rate of 
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transmission, f. 

[0037] For counting sampies of certain amplitude a simple indicator (counting) function defined by equation (1) is 



[0 , x*a 



[0038] The coarse AGO utilizes a coarse gain parameter m c Indicating the amplified output signal strength. For this 
purpose, we Introduce a subset S c of S, representing the set of digital words or quantizer levels that are relevant for 
indicating signal strength. This parameter m e Is obtained by equation (2): 



mc=IX8 ai (bj) (2) 
ajeS 0 H 

bj is the quantized sample at discrete time Index J. For many channels a good choice Is to use S c ={a 2 „}, In contrast, 
for optical channels or generally for binary channels with signal-dependent noise S^a^,,} provides satisfactory 
performance. Note that the latter choice can also be applied for other channels, without loss in performance. 
[0039] During any coarse adjustment phase the coarse logic unit 9 compares the parameter m c with one or more 
thresholds and generates a signal for reducing or increasing the coarse gain accordingly. Consequently, the coarse 
AGC ensures that the signal has sufficient power, in the specified limits. 

[0040] The fine AGC 2 will be adjusted by a uniformity parameter m u , indicating the uniformity of the distribution of 
discrete amplitudes. To this end the set Sof all quantization levels Is partitioned into I disjoint subsets (equation (3)), 
with 2£l£2 n : 



[0041] To each subset S q an integer number k q Is associated, kq represents the expected number of times a quan- 
tization level In the subset S q will be expected In M observations or samples. Since each sample Is associated with 
exactly one subset equation (4) holds: 



Sk q =M (4) 



[0042] Normally, to achieve uniformity, kq should be chosen as W\, With additional knowledge about the channel, 
better performance may possibly be achieved wfth a slightly non-uniform choice of values kq. 
[0043] The uniformity parameter m u , controlling the fine AGC unit, is defined by equation (5): 



-£[xfM>j)U 

q=^B|eS q j=1 j 



(5) 



[0044] With suitable selection of parameters, In particular kq, minimum BER can be achieved by minimizing the 
uniformity parameter m u . An inappropriate selection of parameters can iead lo local minima, which make it hard to find 
the best gain value, producing the lowest number of errors. 
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[0045] Figure 2a shows a possible Implementation of fine logic unit 8. Figure 2a shows a distribution circuit 50, 
counters 51 to 54, absolute difference circuits 69 to 62, registers 55 to 58, adder 70, comparator 72, holding circuits 
71 and 74, exciuslve-OR (XOH) gate 73 and fine gain calculation circuit 75. The digital words provided by ADC 6 on 
line 15 or digital equalizer on line 20 are Input Into distribution circuit 50. Distribution circuit 50 compares each digital 
word to a plurality of digital values the digital words may assume in order to associate each digital word to one of the 
subsets S q , 1 £ q £ 1. To each of the subsets S q one of the counters 51 to 54 Is associated, Distribution circuit 50 
generates a clock pulse for every digital word it receives in order to increment the counter associated with tho subset 
to which the digital word belongs. Additionally, for each subset one absolute difference circuit and one register Is 
provided, Each of the absolute difference circuits calculates the absolute value of the difference of the counter value 
minus the register value of the counter and the register to which the absolute different circuit is connected to. In each 
of the registers one number k^, 1 <, q <l is stored. Adder 70 adds the outputs of all absolute different circuits. Conse- 
quently, distribution circuit 50, counters 51 to 54, absolute difference circuits 59 to 62, registers 55 to 58 and adder 70 
Implement equation (5). Therefore, the output of adder 70 constitutes the uniformity parameter nv 
{0046] The counters 51 to 54 count M digital words input Into distribution circuit 50 during a period of MT wherein T 
Is the symbol period. Then a uniformity parameter mjp] Is calculated at time t p and the counters are reset. 
[0047] Holding circuit 71 stores the previously calculated uniformity parameter m u [p-1] which has been calculated 
at time t^tp-MT wherein t p Is the present time. In addition, the calculation of the uniformity parameter m u may take 
some idle time Tj during which digital words are ignored on lines 15 or 20. If T, Is different from 0,T|must be added to 
MT but, for simplicity, It is assumed to be 0 In the following discussion. Comparator 72 compares the newly calculated 
uniformity parameter mjp] with the previously calculated uniformity parameter mjp-1 J In order to output control signal 
c[p} £{-1,1}. Control signal c[p] is input Into XOR circuit 73. it generates a timing signal bjp] 6 {—1 ,+1} to be used for 
gain phase adjustment. Timing signal b[p] is feed back to the other Input of XOR circuit 73 via holding circuit 74 in order 
to delay timing signal b[p] by MT. Consequently the XOR circuit 73 XORs c[p] and b[p-1 j. Table 1 depicts Its functionality. 



functionality of XOR circuit 73 


c[ P ) 


b[p-1] 


b[p] 




-1 


1 


- 1 


1 


-1 


1 


-1 


-1 


1 


1 


1 



[0048] Practically, minimizing the uniformity parameter m u and consequently the gain adjustment works in a dithering 
mode. During operation this circuit adjusts the gain in a way that the uniformity parameter m u dithers around its mini- 
mum. 

[0049] The new value of the fine gain control signal is obtained by using equation (6): 

fflIp] = f8fp-1]+b[p]A (6) 



[0050] In equation (6) fg[p-1] represents the previous orold fine gain control signal and A represents an adjustment 
step of the fine gain control signal. Equation (6) Is implemented In fine gain calculation circuit 75. The fine gain control 
signal fg is provided via line 1 8 to fine VGA 3 which updates th e gain , After some fine gain adjuslments, the gain range 
wiil be reached, and subsequently tracked, which corresponds to the minimum uniformity parameter m u and BER. 
[0051] In an alternate embodiment XOR circuit 73 may be a conventional logic circuit Inputting either logic 1s or 0s 
and outputtlng Os or 1s. In this embodiment equation 6 must be modified In that previous fine gain control signal fgjp- 
1] Is increased by A if bjp] Is 1 and previous fine gain control signal fg[p-1) is decreased by A if bjp] Is equivalent to 0 
In fine gain calculation circuit 75. 

[0052] Fig. 2b shows a second embodiment of fine logic unit 8. in contrast to the first embodiment shown in Fig. 2a 
the elements 55 !o 75 are replaced by holding circuits 81, 82, 83, 84 and microprocessor 85. Upon processing M 
samples by distribution circuit 50 and counters 61 to 54, the counter values are stored In holding circuits 81 to 84. Then 
the counters are reset and are ready for processing the next M samples. Microprocessor 85 then reads out the values 
of the holding circuits 81 to 84 and calculates a uniformity parameter m u in compliance with equation (5). Then the 
microprocessor performs the functionality of holding circuit 71 , comparator 72, XOR circuit 73, holding circuit 74 and 
fine gain calculation circuit 75 as described in connection with F)g. 2a. 
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[0053] As will be explained later In connection with the simulation result a typical value for M is 120000. As a con- 
sequence the required operating frequency of distribution circuit 50 and counters 51 to 54 is substantially higher than 
the operation frequency of microprocessor 85. Although microprocessor 85 has to perform not only a single operation 
In order to calculate equation (5) and perform the min Imum search a substantial difference in operation speed of roughly 
s afactor 1200 remains for M equivalent to 120000 provided that the processing of the counter values takes 100micro- 
processor clock cycles. So even for high-speed applications a general-purpose microprocessor 85 having a moderate 
clock frequency can be used. The application of standard microprocessors reduces development costs and receiver 
costs especially for small batch sizes. 

[0054] The holding circuits 81 to 84 may be omitted If microprocessor 85 Is quick enough to enter all counter values 
10 within a short time or a certain number of samples are discarded In an IdJe period T, between the blocks of M samples 
which are counted. In the idle period the microprocessor 85 may read out the counter values. 
[0055] In another embodiment the counter values are read out by microprocessor 85 one after another in a distance 
smaller than MT/l. Then each counter value is reset, in this embodiment the periods of time each counter evaluates 
are shifted with respect to each other. Since it is assumed that the symbols do not appear In a certain pattern but rather 
is are statistically Independent this shift may be acceptable. 

[0056] In a third embodiment, which is not shown in a figure, the quantized samples may be directly input into a 
microprocessor which performs the function of elements 50 to 75 shown in Fig. 2a. Such a circuit structure can be 
applied only to moderate transfer rates due to operation speed limitation of Ihe microprocessor 85. In this embodiment 
counters 51 to 64 may be Implemented simply by registers which are incremented by an adder logic. 
so [0057] In a further embodiment subsampling may be performed, i. e. not each sample but rather only e.g. each 
second or third sample Is processed and counted. 

[0058] Figure 3 shows a statistic of digitized sample values in the form of a bar graph. The possible digital words 21 
are plotted on the x-axls. It Is assumed that ADC 6 performs a 3-bit digitization. As a consequence the possible digital 
words range from 0 to 7. The area 23 of each bar illustrates the number of occurrences of a specific digital word wllhln 

25 a sequence of M digital words. Such a statistic could be obtained If one counter Is provided for each digital word, 
[0059] Figure 4 illustrates counter values for the statistics of Figure 3 for a lirst subset. The set S of all possible digital 
words is subdivided into four subsets 22. The lirst subset S, comprises the digital words 0 and 1 , the second subset 
S 2 the digital words 2 and 3, the third subset digital words 4 and 5 and the fourth subset the digital words 6 and 7. The 
area of tho bars having a solid upper boundary represents the counter value of a counter for the respective subset. As 

30 a consequence the area of the bar for subset S, In Figure 4 should be equivalent to the area of the bars for digital 
words 0 and 1 In Figure 3. The broken lines 24, more specifically the area of the bars defined by the broken lines, 
Illustrates the constant values k-,, k 2 , k 3 , k 4 which are subtracted from the counter values In order to obtain difference 
values. As oxplalncd in connoction with equation (5) the absolute vaiues of the difference values are added in order 
to obtain the uniformity parameter m u . The absolute values of the difference values are represented by the area of the 

3s stripes 26 between the solid and the broken lines. 

[0060] Figure 5 Illustrates the choice of two subsets in a further embodiment of this invention. The set S of possible 
digital vaiues is subdivided into an outer subset S,, 41 and an inner subset S z , 42. The outer subset S t comprises the 
digital words 0, 1 , 6 and 7 whereas the Inner subset S 2 comprises the digital words 2 to 5. 11 may be said that the inner 
subset is chosen in a manner that the difference between any digital word comprised In the Inner subset and the 

•to average of all digital words Is smaller than the absolute difference between any digital word comprised in the outer 
subset and the average. The average of all digital words in the example explained in connection with Figures 3 to 5 is 
4, A uniformity parameter in this embodiment may be calculated by equation ('/): 



[0061] As In Figures 3 and 4, the area of bars 43 and 45 represents the count of digital words within the first subset 
and the area of bar 44 represents the count of digital words In the second subset. The broken lines 46 represent an 
offset fy. The uniformity parameter m u defined by equation 7 may assume both, positive and negative values. In a 
preferred embodiment in each adjustment step, the fine gain Is changed by a value proportional to the uniformity 
55 parameter. Provided that the offset value k 1 Is not too big, this is tantamount to saying that the amplification is increased 
if the count associated with the inner subset is greater than the count of the outer subset. 

[0062] The uniformity parameter defined in equation 7 is suitable for being processed by a proportionat-integral- 
differentia! controller or a fuzzy controller since no minimum search Is required. 



45 




(7) 



a; 6 T 2 j=1 
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Simulation 

[0063] The performance of the method has been simuiated for the case of an optical channel distorted by chromatic 
dispersion and noise. In this example, a MLSE (Maximum-likelihood sequence estimator) algorithm has been applied 
for digital equalization, and the parameters utilized in simulation have been selected as: 

M=12O'O0O 

three bit quantization, n=3, 2"=8 
S^fa^ag} 

m^M/8 (M/8=15'000) 

S = S 1 uS 2 uS 3 u S 4 , SHaiM S 2 ={a 3 ,a 4 ), S 3 ={a s ,a 6 } and S 4 ={a 7 ,aa} 
k^M/4 (M/4=30'000), 1 < q < 4 

[0064] The simulation results are shown in Fig. 6, In this graph, the coarse gain parameter m c , the uniformity pa- 
rameter m„ and the BER are plotted against the normalized signal power. In the low signal power regime e. g. beiow 
0.2 normalized signal power, the digital equalizer 7 produces high BER due to existing unoccupied outer quantization 
levels. By Increasing the signal power, the sequence estimator arrives at the region of the best BER for normalized 
signal power between 0.6 and 1.2. For the normalized signal power between 0.6 and 1.2 the parameter m u reaches 
the minimum value and the parameter m 0 is greater than M/8-1 5'uOO. Further Increasing the signal power leads to the 
quantization overloading that causes an Increase of the BER. 

[0065] As explained above, the gain of the coarse VGA unit 2 is adjusted first by increasing the coarse gain from 
lower to higher values. After the m e exceeds M/8 the coarse gain Is kept constant and the fine AGC unit minimizes the 
uniformity parameter m u and provides the best signal power indicating the minimum BER. 
[0066] For unknown but distorted channel characteristics, the existence of the m u minimum enables the AGC to 
deliver useful amplitude statistics for the downstream digital equalizer. Problems of signal-to-noise ratio reduction with 
undistorted, almost ideal channels are obviated by the coarse AGC and especially by the lower bound on m c , In this 
case, the coarse AGC will keep the signal power In some range and the fine AGC will be switched off practically. With 
an undistorted signal detector, e.g. based on amplitude histogram analysis, or by external configuration command, the 
fine AGC could really be switched off with the advantage of less dither penalty. Nevertheless, for time-variant channels 
with distortion the fine AGC demonstrates Its performance, 

[0067] Further modifications and variations of ihe present Invention will be apparent to those skilled In the art In view 
of this description. Accordingly, this description is to be construed as Illustrative only and is forthe purpose of teaching 
those skilled in the art the general manner of carrying out the present Invention. It Is to be understood that Ihe forms 
of the Invention shown and described herein are to be taken as the presently preferred embodiments. 



Claims 

1. A method for controlling amplification within a receiver (1) comprising: 

setting the gain of a variable gain amplifier (3) which amplifies an Input signal (10) to obtain a gain-controlled 
signal (12); 

analog-to-dlgital converting (6) said gain controlled signal (12) thereby generating digital words (15, 20) each 
digital word having a value of a plurality possible digital values; 

counting digital words having a value within a first subset (S^ of possible digital values thereby generating a 
first counter value (51); said counted digital words belonging to a predetermined number of digital words; said 
first subset (S t ) comprising at least one value; 

characterized by 

counting (52, 53, 54) digital words having a value within a second subset (S 2 , S 3 , S 4 ) of possible digital values 
thereby generating a second counter value; 

said counted digital words beionglng to said predetermined number of digital words; said second subset being 
different from said first subset; said second subset comprising at least one value; and 
setting said gain In accordance with the first and second counter value in a fashion that all counter values are as 
equal as possible. 
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The method of claim 1 , wherein each subset comprises at feast two possible digital values. 
The method of claim 1 or 2, said setting further comprising: 

summing all counter values (70) in order to obtain a uniformity parameter; and 

setting (8, 18) said gain in order to minimize said uniformity parameter. 

The method of claim 1 or 2, said setting further comprising: 

subtracting a value (24, 55, 56, 57, 58) from each counter value (25, 51,52, 53, 54) In order to obtain difference 
values (26); and 

summing the absolute values of all said difference values tn order to obtain a uniformity parameter; and 

setting (8, 18) said gain in order to minimize said uniformity parameter. 

The method of claim 1, characterized In that the absolute dilference between any value within said first subset 
(42) and an average of all possible digital values Is greaterthan the absolute difference betv/een any value within 
said second subset (41) and said average and 
said setting said gain further comprising: 

increasing said amplification with the count of said first counter (51) is greater than the count of said second 
counter (52, 53, 54); and 

decreasing said amplification If the count of said first counter is smaller than the count of said second counter, 
The method of claim 5, said sotting said gain further comprising: 

calculating the difference of the first counter value minus the second counter value; 

adding an offset to said difference In order to obtain a uniformity parameter; 

said gain is Increased or decreased by a value proportional to said uniformity parameter. 
A circuit to be used within a receiver comprising: 

a variable gain amplifier (3) for receiving an analog signal and outputtlng a gain controlled signal (12); 

an analog-to-digital converter (6) connected to said variable gain amplifier (3) for digitizing said gain controlled 
signal, thereby generating digital words, each digital word having a value out of a plurality of possible digital 
values; 

a first counter (51 ) connected to said analog-to-dlgital converter for counting digital words having a value within 
a first subset of possible digital values; said counted digital words belonging to a predetermined number of 
digital words; said first subset (S,; 41) comprising at least one value; 

characterized by: 

a second counter (52, 53, 54) connected to said ana!og-to-digita I converter (6) forcountlng digital words having 
a value within a second subset (S 2 , S 3 , S 4 ; 42) of possible digital values; said counted digital words belonging 
to said predetermined number of digital words; said second subset (S 2 , S 3 , S 4 ; 42) being different from said 
first subset (S,; 41); said second subset comprising at least one value; and 

a gain calculation circuit (75) for setting said gal n in accordance with the values of said first and second counters 
in a fashion that said values of said counters are as equal as possible. 

A circuit of claim 7, further comprising an adder (70) which Is connected to each of said counters for adding the 
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counter values; said adder (70) providing its output to said gain calculation circuit (75) which sets said gain in 
accordance with the output signal of said adder (70). 

9. The circuit of claim 7, further comprising: 

a first register (55); 

a second register (56, 57, 58); 

a first absolute difference circuit (59) being connected to said first counter (51) and said first register (55) for 
calculating the absolute difference between the value of said first counter (51) and the value stored in said 
first register (55) 

a second absolute difference circuit (60, 61 , 62) being connected to said second counter (52, 53, 54) and said 
second register (56, 57, 58) for calculating the absolute difference between the value of said second counter 
(52, 53, 54) and the value stored In said second register (56, 57, 58); and 

an adder (70) connected to each of said absolute difference circuits (59, 60, 61 , 62) for adding said absolute 
differences; said adder (70) providing its output to said gain calculating a circuit (75) which sets said gain in 
accordance with the output provided by said adder (70). 

10. A circuit of claim 7, characterized in that the absolute difference between any value within a said first subset (41) 
and an average of all possible digital values smaller than the absolute difference between any value within said 
second subset (42) and said average and said circuit further comprising a subtractorfor calculating the difference 
of said first counter value minus said second counter value; 

said subtractor providing its output to said gain calculating circuit (75) for setting said gain In accordance with the 
■ output of said subtractor. 
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SYMBOL TIMING RECOVERY METHOD FOR LOW 
RESOLUTION MULTIPLE AMPLITUDE SIGNALS 

5 COPYRIGHT NOTICE 

A portion of the disclosure of this patent document contains material which is subject to 
copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone 
of the patent document or the patent disclosure, as it appews in the Patent and Trademark Office 
1 0 patent file or records , but otherwise reserves all copyright rights whatso ever. 

BACKGROUND OF THE INVENTION 

Field of Invention 

15 This invention relates to recovery of symbol tuning of a signal received in a digital 

communication system. The invention is more particularly related to the detection of symbol edges 
for use in symbol timing recovery of low resolution digital signals encoded in an analog received 
signal. 

20 Discus sion of Back ground 

hi general digital communication systems, a transmitting apparatus modulates a digital signal 
(1) by any known modulation techniques (Qiiadrauire Amplitude Modulation (QAM), or Quadrature 
Phase Shift Keying (QPSK), for example). The modulation technique conveits the digital signal into 
an analog signal for transmission. A receiver of the analog modulated signal demodulates the signal 

25 and converts it into a digital signal, restoring it to the original form (digital signal (1)). 

Many modulation and demodulation tecliniques ate biow that are capable of accurately 
restoring the received signals into original digital signals (also referred to as symbol values) when 
the received signals are sampled at the same interval of time (frequency) and symbol starting point 
(phase) as those of the symbols transmitted from the transmitter. Without proper frequency and 

30 phase alignment in the receiver (Rx), increased computations are required to maintain accuracy of 
the digital signal restoration, or additional errors are encountered (increasing a bit error rate (BBR) 
of the communications system in general, and the receiver in particular). Therefore, one of the 
important processes that occur in the receiver of the digital communication system is to accurately 
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restore a symbol timing (frequency and phase) representing a point in time that the received signal 
should be sampled. 

In some systems, a specific carrier corresponding to a sampling frequency is transmitted, 
providing freqxiency and phase for the sampling. However, many modem systems do not provide 
5 a carrier, and it is therefore important to accurately recover a symbol timing using only the received 
signal, "With the proliferation of digital communication systems, there has been an increased need 
for improved and distinctive symbol timing recover (STR) processes within receivers of digital 
communications. A conventional apparatus for recovering a symbol timing is shown in FIG. 1 . 

In FIG. .1, a received signal is input to an analog-to-digital (A/D) converter 100. The A/D 
10 converter 100 samples the input signal at a frequency supplied from a voltage controlled oscillator 
(VCO) 110 and converts the analog signal into a digital signal. The digitized signal is output via a 
matching filter 120, and simultaneously input to a timing error detector 130. The timing error 
detector 130 detects a symbol timing error from the input signal and outputs a detected error value, 
to a loop filter 140, The loop filter 140 filters the input signal to remove various noises, The VCO 
15 110 oscillates at a frequency according to the output of the loop filler 140 and supplies the oscillated 
frequency signal to the A/D converter 100. 

However, in conventional systems, receiver (Rx) and transmitter (Tx) inaccuracies, operating 
conditions, and other reasons can cause symbol timing drift that results hi the number of symbols 
transmitted from the transmitter differs or has significantly different tuning from that of the restored 
20 signals. For example, one reason is because the receiver does not know when the symbols will 
anive, which can be due to a number of factors, including propagation delay between the Tx and Rx, 
particularly in a wireless cormniuiication environment where distances may become large. Another 
reason is that the transmitter and receiver clocks drift from their nominal values, forcing the receiver 
to continually adjust (or adapt) to maintain the appropriate sample time instance. By appropriate, 
25 it is desired that the sampling time instance match (or closely match) a maximum eye opening of the 
received signal, since this timing epoch corresponds to a maximum signal to noise ratio (SNR) value 
of the received signal. 

One common method for continually adjusting to maintain the appropriate sample time 
instance is known as matched filter technique, which is illustrated in Fig. 2. In Fig. 2, a fixed 
3 0 oscillator 200 generates a constant frequency signal and supplies the same to an A/D converter 210. 
The A/D converter 210 converts the received signal into a digital signal using the frequency signal 
supplied from the fixed oscillator 200, and outputs tire digitized signal to an interpolator 220. 
Meanwhile, a controller 230 receives a symbol timing error signal output from a timing error 
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detector 250 via a loop filter 240, The controller 230 generates an interpolation coefficient according 
to the input symbol timing error signal and supplies the same to the interpolator 220. The 
interpolator 220 interpolates the digitized signal with the interpolation coefficient supplied from the 
controller 230, and outputs the interpolated result to a matching filter 260. 
5 The A/D converter 210 samples the received signal at a sampling frequency supplied from 

the fixed oscillator 200 and converts the analog signal to the digital signal. The digitized signal is 
output via the interpolator 220 and the matching filter 260 and simultaneously input to the timing 
error detector 250. The symbol timing error signal which is detected in the timing error detector 250 
and has noise removed in the loop filter 240 is then input to the controller 230. The controller 230 

10 generates the interpolation coefficient according to the input signal. The interpolator 220 
interpolates the digitized signal according to the interpolation coefficient supplied from the controller 
230, and thereby removes the symbol tuning error. Although this apparatus can accurately recover 
the symbol timing, it is very complicated to implement using hardware. 

However, none of these techniques provide the best solution for modern digital 

15 communications devices. Moreover, the above techniques favor high resolution received signals, 
In this application, resolution refers to discrete time representations of a continuous signal. High 
resolution refers to a relatively large amount of bits in the ADC. For example, a number of bits (b), 
such as b >= 6, are considered high resolution. 

20 SUMMARY OF THE INVENTION 

The present inventors have realized that modem digital commiuiications systems fail to 
address the needs of applications of low resolution digital systems. 

The present invention provides a robust method with good BER performance for low 
resolution signals. The present invention can be implemented without the complexities seen in the 
25 ' prior art, and can be modified for improved performance when high resolution signals are available. 

Roughly described, the present invention provides Symbol Timing Recovery (STR) via a set 
of edge detections maintained in a histogram that are utilized to calculate the phase and timing of 
subsequently received symbols. 

The present invention is embodied as a device, comprising, an ADC mechanism configured 
30 to sample a received signal at a predetermined over-sampling rate, a detector connected to said AD C 
mechanism and configured to detect a symbol identified by at least one of said samples during 
predetermined symbol times, and a symbol timing mechanism connected to said ADC mechanism 
and said detector and configured to read said samples and adjust a timing of the detection performed 
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by said detector based on previously read samples so that the detected samples are at a most likely 
maximum eye opening of said samples during said predetermined symbol times. 

The present invention includes a method of symbol timing, or recovery of symbol timing, 
comprising the steps of, sampling a received signal at a predetermined over-sampling rate, 
5 detemiining a timing for detection of a specific one or set of said samples based on historical data 
of samples previous samples, and detecting said specific one or predetermined number of samples 
based on the determined timing. 

Both the device and method may be conveniently implemented on a general purpose 
computer, or networked computers, and the results may be displayed on an output device connected 
10 to any of the general purpose, networked computers, or transmitted to a remote device for output or 
display. 

ftRJEF DESCRIPTION OF THE DRAWINGS 

A more complete appreciation of the invention and many of the attendant advantages thereof 
15 will be readily obtained as the same becomes better understood by reference to the following detailed 
description when considered in connection with the accompanying drawings, wherein: 

Fig. 1 is conventional apparatus for recovering a symbol timing; 

Fig, 2 is a block diagram illustrating a conventional matched filter technique for continually 
adjusting to maintain an appropriate sample time instance; 
20 Fig. 3 is a block diagram illustrating the Symbol Timing Recover section in a digital receiver; 

Fig. 4A is a graph illustrating a waveform of a digitally sampled continuous analog signal; 

Fig. 4B is a graph illustrating a waveform of a digitally sampled low resolution continuous 
analog signal; 

Fig, 5 is a flow chart illustrating a detector adjustment process according to the present 
25 invention; 

Fig. 6 is a block diagram of a symbol edge detection device according to one embodiment 
of the present invention; 

Fig. 7 is a flow chart illustrating one embodiment of a process for utilizing histogram data 
for determining a most likely eye opening for a next symbol time; 
30 Fig. 8 is a flow chart illustrating one embodiment of mamtaining a ninning average of timing 

for a most likely max eye opening of a next symbol time; 

Fig. 9 is a block diagram of a symbol center value detection device according to another 
embodiment of the present invention; 
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Fig. 10 is an example timing chart of symbol and e value determinations according to the 
present invention; 

Fig. 11A is a set of histogram plots for varying numbers of symbols (M) developed for an 
Rx eye diagram; 

5 Fig. 11B is a 4 level eye diagram, prior to ADC, using a flex-paging protocol for 

demonstration purposes; 

Fig. 12 is a set of receiver Eye diagrams and corresponding histograms for analyzing 
performance of a rising edge embodiment of the present invention with a low SNR signal; 

Fig, 13 is a set of receiver Eye diagrams and corresponding histograms for analyzing 
1 0 performance of a center detection embodiment of the present invention with a low SNR signal; 

Fig. 14 is a set of graphs illustrating BER vs. Bb/No (dB) performance of various 
embodiments of the present invention vs. ideal and theory STR's; 

Fig. 15 is a set of graphs illustrating BER vs, Eb/No (dB) performance of various 
embodiments of the present invention vs. ideal and theory STR's, including 2-FM (binary frequency 
15 shift keying) and 4-FM (4 level frequency shift keying) theory. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In a digital communications system, a number of bits used in the Analog to Digital Converter 
(ADC) comprises a tradeoff between power consumption, board real estate, and performance, In the 
20 present invention, low resolution signals are targeted (for example, signals used in pagers). An 
example low-resolution system would, for example, choose b = 2 bits for the ADC and 4-level 
baseband signals (low frequency content signals with discrete amplitude levels and carrier removed, 
for example). 

Referring again to the drawings, wherein like reference numerals designate identical or 
25 corresponding parts, and more particularly to Fig, 3 thereof, there is illustrated a block diagram 
illustrating the Symbol Timing Recover section in a digital receiver. A received signal is fed to an 
ADC 300 which is timed via an oscillator 310. The digitally converted signals b-bits 315 are fed to 
a timing recovery module 320 and a detector 330. The timing recovery module 320 identifies the 
symbol edge and adjusts the phase of detector 330 such that the detector 330 recognizes the symbol 
30 at an amplitude sufficient to recover a bit value previously represented in the received signal. The 
recovered bits, or symbols, are output. 

Fig, 4A is a graph illustrating a waveform of a digitally sampled continuous analog signal 
(for 8x over-sampling), hi this example, over-sampling is the number of digital bits converted from 
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the received signal for a length of time equivalent to a period of a single bit in the digital waveform 
to be recovered. Thus, fig. 4A shows 8 samples, at times y 0j yi, ya> yj, y*, y 5 , y<s, and y 7l each 
provided by the ADC, from which a symbol is to be recovered. For the eye diagram of fig; 4A, the 
detector would want to detect the symbol at t - y% because V3 represents the time at which the eye 
5 is at maximum (greatest amplitude of the symbol to be recovered). 

When using a matched filter approach, each digitally sampled waveform (Fig. 4A, or 4B, for 
example) is compared to a stored waveform that represents (or "looks" like) a waveform that is 
expected to be received. Each sampled waveform is compared, and a stored waveform having the 
highest correlation to the sampled waveform is chosen for output 

10 Fig. 4B is a graph illustrating a waveform of a digitally sampled low resolution continuous 

analog signal (for 8x over-sampling). The digitally sampled signal shows a flattening effect 420. 
The lost information (additional amplitude at the sample points missing due to low resolution) 
makes matching more difficult because less variance between signals makes each signal have less 
correlation to the stored waveforms. 

1 5 Therefore, it can be seen that the matched filter technique works well with high resolution 

signals. However, with low resolution signals, performance is degraded because much of the 
information contained in the signal is lost due to resolution reduction. 

The present invention observes symbol (or bit) edges in the received digitally sampled signal, 
and develops a histogram (normalized histogram, or probability density function, normalized by the 

20 total number of events, for example) of the boundaries of each symbol. Fig. 5 is a flow diagram that 
shows the basic process. At step 500, a sample is taken (MSB, for example). At step 510 a second 
MSB sample is taken. The samples (MSB's in this example) are compared to determine if a 
transition has occurred (step 520). A transition is detected by determining if Sample 1 (SI) and 
Sample 2 (S2) are different (different samples indicates a transition, equivalent, or roughly 

25 equivalent samples indicate no transition). 

The detection information is used to develop a histogram that identifies timing and recovery 
of past edges or transitions of the sampled received signal, With each transition, the histogram is 
updated (step 530). The histogram is utilized to adjust the detector (step 540). The detector can be 
adjusted to any of edge, center, or other portions (1/4, for example) of the incoming symbols, 

30 depending on the type of detector used, To finish the loop, SI is substituted with SI and the loop 
continues (step 550). Variations of the flow illustrated in Fig. 5 can be made based on the present 
disclosure. 
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Fig, 6 illustrates a block diagram of a symbol edge detection device according to pne 
embodiment of the present invention, m Fig. 6, ADC 600 provides an MSB and a LSB recovered 
by sampling of a received signal (b «= 2). 

An edge detector 610 compares adjacent MSB samples, The edge detector 610, in this 
5 example, is constructed from a delay circuit 612 and an exclusive or circuit device 614. 
Alternatively, the edge detector may be constructed from other configurations of parts, including, 
for example, a delay circuit and a comparator. What is important is that the edge detector is capable 
of determining a transition of the MSB (or other sampled bit), The edge detected (e) is then fed to 
an N-bin histogram 620. In one embodiment, the n-bin histogram 620 is constructed of a memory 

10 device that maintains a history of timing for each edge detected. In another embodiment, one or 
more averages of edge detectors are calculated and stored. 

A detector 640 detects a bit in the sampled received signal at a point that is most likely to be 
the highest amplitude of the bit being recovered from the received signal. A center adjust device 630 
utilizes the histogram 620 to determine the point where the highest amplitude occurs based on 

15 previous edge detections. For example, if the last n edge detections each came in at a same time 
instance interval, that same interval would be the most likely next edge and the most likely point of 
highest amplitude would be calculated based on that interval. Using specific values, if each symbol 
is received in .25 ms intervals, the next highest amplitude point could bo calculated at ,25ms + Ts/2, 
where Ts = the timing period of the received bits, 

20 In other words, the edge detector compares adjacent MSB samples. If the adjacent samples 

are the same, then e = 0; which indicates that a transition has not occurred. If the adjacent samples 
are different, then e=4; indicating that an edge has been detected, When an edge is found, the time 
instances (time samples from where the edge was found) are distributed into the N-bb histogram 
620. In one embodiment, N=8 represents the over-sampling ratio of the ADC (the number of 

25 samples per bit or symbol). 

Thus, statistics for an edge detection are gathered. The statistics are accumulated for a pre- 
detemiined number of Symbols M(M-40, for example). Using the accumulated statistics, a most 
likely (ML) position of a symbol edge can be determined. For example, given a set of accumulated 
statistics, and assuming an eye that is at least somewhat symmetrical, the center adjustment is made 

30 such that eye is assumed to be at a detected edge, plus Ts/2, Ts being the symbol time (period of the 
symbols being decoded). Ts is calculated based on the histogram data. Ts is calculated using any 
method, peak-to-peak, for example. Ts, generally being constant, is known ahead of time. However, 
it is not known where in the histogram the peaks or Ts will fall. In addition, other factors affect 
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timing and placement of Ts, such as propagation time delay, frequency drift of transmitter/receiver, 
interference, etc. Ts may equal T, or be some multiple or fraction thereof (.9T, UT, for example). 

The histogram may be constructed and maintained in many different ways. What is 
important is that a timing of a most likely max eye opening either be maintained or be derivable from 
5 data stored in the histogram. For example, Fig. 7 provides a flow chart of one embodiment of a 
process for using histogram data storing the timing of a max eye opening for each set of samples for 
each symbol time. At step 700, the samples from a sample period are read. The read samples are 
compared to determine the max eye opening (highest amplitude, for example) (step 710). A timing 
of the max eye opening from an edge detection of the symbol time to the max eye opining is stored 

10 in the histogram (step 720) (max eye opening). If all the stored samples are not yet read, samples 
from another period (symbol time) are read and the process repeats (step 740). If all the stored 
samples are read, an average or other method is utilized to determine a most likely max eye opening 
for a next symbol time (step 750). The most likely max eye opening for the next symbol time is then 
added to the tuning of an edge detection of the next symbol time to provide a tuning for the detector 

15 to detect the symbol value. 

Other histogram configurations and processes may be utilized to determine the detection 
timing. Fig. 8 Illustrates another embodiment where a running average of max eye opening timing 
relative to edge detection is maintained. At step 800 sample data is collected for a symbol time of 
a received signal. The collected sample data is analyzed to determine a timing of the max eye 

20 opening of the collected sample data symbol timo (step 810). A running average or other statistic 
representing a most likely max eye opening for a next symbol time is then adjusted based on the max 
eye opening timing of the collected sample data (step 820). The running average is then applied to 
the edge detection of the next symbol tune to provide a detection timing for the symbol. 

Fig. 9 illustrates a block diagram of another embodiment of the present invention. In this 

25 embodiment, a symbol center value is determined. ADC 600 provides an MSB and a LSB recovered 
by sampling of a received signal (b = 2). A transition detector 910 compares adjacent MSB samples. 
The transition detector 910, in this embodiment, is constructed from a delay circuit 912 and an 
exclusive or circuit device 914, Alternatively, the transition detector may be constructed from other 
configurations of parts, including, for example, a delay circuit and a comparator. What is important 

30 is that the edge detector is capable of deterrmning a transition of the MSB (or other sampled bit) (an 
edge triggered device, for example, outputting a fixed value when the input goes from zero to high, 
and outputting zero when the input falls). 
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The value output from the transition detector 910 is e, pulse synched with the symbol 
boundaries. An illustration is provided by Fig. 10. An MSB is illustrated cycling high for a time 
period t| , low for a time period t 2 , and back high for time period t 3 . e is shown pulsing high at each 
transition (1000, 1005, 1010, and 1015). 
5 The edge detected (e) is then inverted by inverter 915. The inverted edge (e bar 1 020) is a 

signal that pulses high at each symbol center. The e bar 1020 signal is fed to a rising edge true 
device 925. The rising edge true device, which operates as a dirac delta type function (impulse type 
function), or one cycle of a fast (fastest) clock cycle, for example, to produce e R 1030. The verified 
rising edge value eR is fed to an N-bin histogram 920, where it is saved. In this embodiment, the 

1 0 center value is adjusted directly, rather than using the edge detection plus l A a symbol time. 

Fig. 11A is a set of histogram plots for varying numbers of symbols (M) developed for a 
received signal (see Rx eye diagram). The SNR for the illustrated diagrams is 1 0 dB. The histogram 
plots provide an analysis of performance of various embodiments of the present invention. 
Histogram for each of M = 10, 20, 40, and 80 are shown. As can be seen from the diagrams, in 

15 general, the edge detections are consistent and provide a reasonable estimate of symbol time. 
However, the eye diagrams show corresponding better averages with increasing numbers of 
histograms averaged (note more distinct peaks, at symbol times 1 and 2, for example, and the 
flattening of non symbol time areas, at approximately .25 and ,75, for example, as the number of 
symbols increases). 

20 Fig. 11B is a 4 level eye diagram, prior to ADC, using a flex-paging protocol for 

demonstration purposes. The present inventors obtained good results using a receiver bandpass filter 
having 8kHz bandwidth, and a post detection fdter having a 2 kHz bandwidth. 

Fig. 12 is a receiver Bye diagram and corresponding histogram for analyzing performance 
of a rising edge embodiment of the present invention with edge detection and a low SNR signal. The 

25 Rx eye chart shows the received signal, and the M=80 histogram illustrates consistent, symmetrical, 
symbol times. 

Fig. 13 is a receiver Eye diagram and corresponding histogram for analyzing performance 
of a center detection embodiment of the present invention with a low SNR signal. Again, the Rx eye 
chart shows the received signal, and the histogram illustrates consistent, symmetrical, symbol times. 
30 In most cases, system performance is equivalent whether using an edge or center detection 

embodiment. However, the system has more jitter kit, then the edge detector histogram will be less 
peaky, mote flat, hi this case, a system designer would want to go with FIG. 13, center detection, 
because jitter means there is uncertainty, if you don't have any jitter, then the designer would want 
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to go with the edge detector, In one embodiment, both center and edge detection is provided, giving 
the option of choosing one versus the other, depending on what the overall received signal is. 

Fig. 14 is a set of graphs illustrating BER vs. EbMo (dB) performance of various 
embodiments of the present invention vs. ideal and theory STR's. As can be seen in the graph, BER 
5 performance for each of N=20, 40, 80, and 1 60 are grouped in close proximity of an ideal SIR. The 
ideal STR is seen grouped with the N= 200..160 group, and deviated from a Theory SIR by 
approximately 1 dB at each portion of the graph. In an implementation of the invention, the number 
of histograms considered can be varied depending on any factor (dB required, or available processing 
power, for example). When additional processing power is available or additional dB needed, M 

10 may be set to ahigher value (160, for example), but in low dB situations (pilot signal, for example), 
a lower value such as 20 may be sufficient. 

Fig. 15 is a set of graphs illustrating BKR vs. Eb/No (dB) performance of various 
embodiments of the present invention vs. ideal and theory STR's, including 2-FMand 4-FM theory. 
Portions of the present invention may be conveniently implemented using a conventional 

1 5 general purpose or a specialized digital computer or microprocessor programmed according to the 
teachings of the present disclosure, as will be apparent to those skilled in the computer art, 

Appropriate software coding can readily he prepared by skilled programmers based on the 
teachings of the present disclosure, as will be apparent to those skilled in the software art. The 
invention may also be implemented by the preparation of application specific integrated circuits or 

20 by interconnecting an appropriate network of conventional component circuits, as will be readily 
apparent to those skilled in the art. 

The present invention includes a computer program product which is a storage medium 
(media) having instructions stored thereon/in which can be used to control, or cause, a computer to 
perform any of the processes of the present invention. The storage medium can include, but is not 

25 limited to, any type of disk including floppy disks, mini disks (MD's), optical discs, DVD, CD- 
ROMs, micro-drive, and magneto-optical disks, ROMs, RAMs, EPROMs, EEPROMs, DRAMs, 
VRAMs, flash memory devices (including flash cards), magnetic or optical cards, nanosystems 
(including molecular memory ICs), RAID devices, remote data storage/archiveAvarehousing, or any 
type of media or device suitable for storing instructions and/or data. 

30 Stored on any one of the computer readable medium (media), the present invention includes 

software for controlling both the hardware of the general purpose/specialized computer or 
microprocessor, and for enabling the computer or microprocessor to interact with a human user or 
other mechanism utilizing the results of the present invention . Such software may include, but is not 
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limited to, device drivers, operating systems, and user applications. Ultimately, such computer 
readable media further includes software for performing the present invention, as described above. 

Included in the programming (software) of the general/specialized computer or 
microprocessor are software modules for implementing the teachings of the present invention, 
5 including, but not limited to, comparing consecutive MSB values, maintaining a iiistogram of signal 
samples, averaging max eye openings of sample sets, adjusting detection times, and the display, 
storage, or communication of results according to the processes of the present invention. 

Obviously, numero\is modifications and variations of the present invention are possible in 
light of the above teachings. It is therefore to be understood that within the scope of the appended 
1 0 claims, the invention may be practiced otherwise than as specifically described herein. 
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WHAT IS CLAIMED IS : 

1 . A symbol timing device, comprising; 

an ADC mechanism configured to sample a received signal at a predetermined over-sampling 

5 rate; 

a detector connected to said ADC mechanism and configured to detect a symbol identified 
by at least one of said samples during predetermined symbol times; and 

a symbol timing mechanism connected to said ADC mechanism and said detector and 
configured to read said samples and adjust a timing of the detection performed by said detector based 
10 on previously read samples so that the detected samples are at a most likely maximum eye opening 
of said samples during said predetermined symbol times. 

2. The symbol timing device according to Claim 1 , wherein : 
said symbol timing mechanism comprises, 

15 a histogram device configured to, 

store a histogram of said samples over a predetermined number of said symbol times, and 
calculate a most likely max eye opening of said samples of a next symbol time based on said 
histogram, and 

an adjustment device configured to adjust said liming of the detection performed by said 
20 detector based on the calculated most likely max eye opening. 

3 . The symbol timing device according to Claim 2, wherein: 
said over-sampling rate is N; and 

said histogram device comprises, 
25 a memory configured to store M sets of samples, each sample set comprising N samples, and 

each sample set corresponding to one of said predetermined symbol times, and 

a calculation mechanism configured to detennine said most likely eye-opening of said set of 
samples of a next symbol time based on the sample sets stored in said histogram. 

30 4, The symbol timing device according to Claim 2, wherein said histogram device 

comprises a calculator configured to calculate said most likely max eye opening by averaging max 
eye openings of said samples over said predetermined number of symbol times. 
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5 . The symbol timing device according to Claim 1 , wherein : 
said symbol timing mechanism comprises, 

an edge detector configured to deteimine an edge of symbols sampled by said ADC 
mechanism; 

a histogram configured to save samples from said ADC for a predetermined number of 
previous symbol times; and 

a center adjustment mechanism configured to adjust a sample detected by said detector based 
on said histogram and edges detected by said edge detector. 

6. The symbol timing device according to Claim 5, wherein : 
said edge detector comprises, 

a delay mechanism having an input connected to said ADC mechanism, and an output, and 
an exclusive or device having one input connected to said ADC mechanism, a second input 

connected to the output of said delay mechanism, and an output connected to said histogram; and 
said histogram uses the output of said exclusive or to identify samples maintained in said 

histogram, 

7. The symbol timing device according to Claim 1 , wherein : 
said symbol timing mechanism comprises, 

a transition detector configured to determine an edge of symbols sampled by said ADC 
mechanism and output a transition detected signal, 

a conversion mechanism configured to convert the transition detected signal to a pulse 
synched with a center of symbols corresponding to detected transitions, and 

a histogram that maintains data on a predetermined munber M of said pulses; 

said detector is further configured to detect symbols based on the data maintained in said 
histogram, 

8. The symbol timing device according to Claim 7, wherein said center adjustment 
mechanism is further configured to adjust time lhat said symbol is detected by said detector.is based 
on an average of past pulses. 
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9. The symbol timing device according to Claim 7, wherein said center adjustment 
mechanism is further configured to adjust time that said symbol is detected by said detector is based 
on a weighted average of past pulses, 

1 0. The symbol timing device according to Claim 7, wherein : 
said conversion mechanism comprises, 

an inverter connected to said transition detector and configured to invert the transition detect 
signal such that the inverted signal is now synched with a center of a symbol that caused the 
transition detection and 

a rising edge detect mechanism connected to said inverter and configured to conveit the 
synched signal to a pulse by said rising edge detect mechanism. 

1 1 . A method of symbol timing recovery, comprising the steps of: 
sampling a received signal at a predetermined over-sampling rate; 

determining a timing for detection of a specific one or set of said samples based on historical 
data of samples previous samples; and 

detecting said specific one or predetermined number of samples based on the determined 

timing. 

12. The method according to Claim 1 1, further comprising the step of: 
adjusting said timing based on newly acquired samples, 

13. The method according to Claim 1 1, further comprising the step of: 

maintaining a histogram of sample sets, each set comprising a predetermined number of 

samples for each of a predetermined number of symbol times; 

wherein said step of determining a timing comprises the steps of, 

averaging a max eye opening of each set of samples maintained in said histogram, and 

calculating a most likely max eye opening for a next set of samples based on the average max 

eye opening. 

14. The method according to Claim 13, wherein said step of averaging comprises performing 
a weighted average of said max eye openings of each set of samples maintained in said histogram. 
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15. The method according to Claim 1 1, further comprising the step of: 

maintaining a histogram comprising a raining average of timing of a max eye opening of 
samples for each of a predetermined number of symbol times; 

wherein said step of determining a timing for detection of a symbol comprises using the 
running average. 

16. The method according to Claim 1 1, further comprising the steps of: 

maiiitaining a histogram indicating timings from an edge detection to max eye opening of the 
sampled signal during past symbol times; 

averaging the timings of past symbol time max eye openings; and 

detecting an edge of a symbol contained in the current sampled signal; 

wherein said step of determining a timing comprises, applying the averaged past symbol max 
eye opening timing to a timing of the detected symbol edge of the current sampled signal. 

17. The method according to Claim 16, wherein said step of detecting an edge comprises the 
steps of: 

delaying the received signal by one symbol time; 

comparing a currently received signal and the delayed signal; and 

signaling an edge detection if die compared signals are at different logic levels. 

18. The method according to Claim 16, further comprising the steps of: 

inverting the detected edge so that the rising edge of the inverted signal corresponds to a 
center of a symbol time of the symbol for which die edge was detected; 
converting the inverted signal into a pulse; and 
wherein: 

said step of determining a timing comprises applying each pulse as a timing for detecting said 
specific one or set of said samples. 

19. The method according to Claim 1 8, further comprising the step of: 
adjusting a timing of said pulse based on data contained in said histogram. 
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20. A symbol timing device, comprising: 

ADC means for sampling a received signal at a predetermined over-sampling rate; 

detector means for detecting a symbol identified by at least one of said samples during a 
predetermined symbol time; and 
5 symbol tuning means for reading said samples and adjust a timing of the detection performed 

by said detector based on previously read samples so that the detected samples are at a most likely 
maximum eye opening of said samples. 

2 1 . The symbol timing device according to Claim 20, wherein : 
10 said symbol timing means comprises, 

histogram means for storing a histogram of said samples and determine a most likely max 
eye opening of said samples based on said histogram, and 

adjustment means for adjusting said timing of the detection performed by said detector 
based on said histogram determined most likely max eye opening. 
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